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We present a large class of three-dimensional spin models that possess topological order with sta- 
bility against local perturbations, but are beyond description of topological quantum field theory. 
Conventional topological spin liquids, on a formal level, may be viewed as condensation of string-like 
extended objects with discrete gauge symmetries, being at fixed points with continuous scale symme- 
tries. In contrast, ground states of fractal spin liquids are condensation of highly- fiuctuating fractal 



objects with certain algebraic symmetries, corresponding to limit cycles under real-space renormal- 
T-H ization group transformations which naturally arise from discrete scale symmetries of underlying 

fractal geometries. A theoretical framework for unified analysis on physical and coding properties 
of fractal spin liquids is presented, which is applicable to all the previously known examples. 
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I. INTRODUCTION 

In a quantum many-body system at zero temperature, topological order may arise when a gapped ground state 
possesses long-range entanglement that cannot be detected by any local measurement or local order parameter [T]. 
Ground state properties of topologically ordered systems are stable against any types of local perturbations, regardless 
of symmetries of perturbations, and depend only on global properties of geometric manifolds on which the whole 
systems are supported [2] ■ The discovery of topological order, such as fractional quantum hall systems [31 [4] , came 
as a great surprise as it was beyond description of the Landau-Ginzburg theory which was once believed to be the 
ultimate theory of a quantum many-body system (see [5] for a review). It is now widely believed that the notion 
of topological order is essential in understanding the emergence of quantum phases with no local order in gapped 
quantum spin liquids, as seen in some resonating valence-bond states of frustrated anti-ferromagnets [6Ul3). In 
addition to fundamental importance in condensed matter physics, the study of topologically ordered systems is of 
practical importance since they are physically natural platforms for realizations of fault-tolerant quantum information 
processing Searches for novel topological phases and their classification in a variety of theoretical and experimental 
systems have been an active area of research. 

For topologically ordered spin liquids with discrete gauge symmetries, their low energy behavior is relatively well 
understood on a formal level as they can be effectively described by topological quantum field theory (TQFT) [21 [7] , a 
field theory with an invariance under continuous deformations (diffeomorphism) [21 [TS] . This is due to an observation 
that physical properties of topologically ordered systems do not depend on local structures of the systems, and 
topological properties of geometric manifolds play a crucial role in classifying topological phases. Indeed, a fairly 
complete class of two-dimensional TQFT-based spin systems with non-chiral topological order has been proposed by 
Levin and Wen where condensation of highly-fluctuating extended objects, called "string-nets" , are responsible for 
emergence of topological order ^16 . A different, but closely related approach based on renormalization group (RG) 
ansatz of ground state wave-functions at fixed points has further verified the validity of TQFT as effective theory of 
topological phases in the presence of continuous scale invariance |17) . 

Yet, in some cases, quantum spin liquids may exhibit topological order that is beyond description of TQFT. For 
example, in three spatial dimensions, the Cubic code, recently proposed by Haah [18,, possesses topological order with 
the stability against local perturbations, but seems to exhibit physical properties that are completely different from 
those of conventional topological spin liquids. For one thing, the number of degenerate ground states is exponential 
in the linear length of the lattice. Furthermore, unlike string-net condensates, the model is free from string-like 
extended objects, and the mobility of quasi-particle excitations is highly constrained via some algebraic rules. It has 
been argued that these exotic features are particularly useful for storing quantum bits (qubits) securely in thermal 
environment, but with some subtlety pointed out due to its thermal instability at non-zero temperature The 
discovery of the Cubic code and relevant models [201 HI] clearly indicates that classification of topological phases via 
TQFT is incomplete; TQFT is just a subset of some universal theory of topological order which is yet to be found. 
The necessary first step is to find a family of topological spin liquids that are beyond TQFT, and to develop a general 
theoretical framework to discuss their physical properties in a unified way. 

The goal of this paper is to present a large class of exactly solvable topological spin liquids on a three-dimensional lat- 
tice which possess exotic topological order that cannot be described by TQFT. Instead of string- like (one-dimensional) 
or membrane-like (two-dimensional) objects with continuous geometries, ground states of the model are condensation 
of extended objects with non-integer dimensionality, namely fractal objects. In this paper, we discuss physical proper- 
ties of such quantum fractal liquids, including; (a) the number of ground states and its dependence on geometric shapes 
of the manifolds, (b) the presence of topological order that is beyond TQFT, (c) propagations of quasi-particle exci- 
tations and their algebraic properties, and (d) topological entanglement entropies as an order parameter. Theoretical 
tools for analyses on quantum fractal liquids are also developed. 

Emergence of fractal objects in correlated spin systems is not a completely new idea. In [22], Newman and Moore 
proposed a toy model of translation symmetric classical spin liquid on a two-dimensional square lattice. It has 
degenerate ground states whose spin configurations resemble the Sierpinski triangle with the famous "log 3/ log 2" 
fractal dimension. Further studies revealed that a family of such fractal systems, refereed to as classical fractal 
liquids, are particularly useful for storing a large amount of classical information reliably as they saturate a theoretical 
limit on classical information storage capacity of local Hamiltonians with mass gap fIT\. Ground states of classical 
fractal liquids do not have continuous scale symmetries, but have discrete scale symmetries only, exhibiting limit 
cycle behaviors under real-space RG transformations as demonstrated in this paper. Physical properties of classical 
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fractal liquids are intrinsically akin to those of a black hole since inner states are completely determined by degrees 
of freedom on a one-dimensional surface, and their coding properties obey an "area law" asymptotically [23j . Such 
exotic features of classical fractal liquids indicate a possibility of novel quantum phases that cannot be described by 
a field theory with continuous scale invariance. Construction of quantum fractal liquids can be viewed as natural 
generalization of classical fractal liquids to a quantum setting. 

We first present a physical picture of quantum fractal liquids in section |ll] by reviewing how condensation of 
extended objects emerge in topological spin liquids. In conventional topological spin liquids described by TQFT, 
extended objects with continuous geometries and integer dimensionality may emerge from underlying discrete gauge 
symmetries [71|5j. In contrast, quantum fractal liquids are condensation of fractal objects with discrete geometries 
and non-integer dimensionality which emerge from certain algebraic symmetries. Geometric properties of extended 
objects can be characterized by topological classes of non-trivial symmetry operators where fractal symmetry operators 
are associated with quantum fractal liquids. In particular, geometric shapes of symmetry operators can be used as 
"order parameters" for classifying quantum phases where discontinuity in quantum phase transition is triggered by 
non-analytic changes of topologically distinct logical operators. 

We then present a general framework to construct a family of classical fractal liquids with a variety of fractal ge- 
ometries in section [nil Quasi-particle excitations in classical fractal liquids may violate charge conservation since they 
do not possess gauge symmetries, but obey more general algebraic symmetries which are associated with polynomials 
that generate fractal geometries. Classical fractal liquids differ from conventional classical spin liquids as correlation 
functions have peculiar oscillatory behaviors associated with power-law decay with imaginary scaling dimensions. We 
demonstrate that ground states of classical fractal liquids correspond to limit cycles of real-space RG transformations. 
In section |IV[ we present a general framework to construct a family of three-dimensional quantum fractal liquids 
by weaving a pair of two-dimensional classical fractal liquids coherently such that their two-dimensional holographic 
images are two-dimensional classical fractal liquids. It is shown that quantum fractal liquids are topologically ordered 
with stability against local perturbations, but are beyond TQFT. We then discuss algebraic properties of quasi-particle 
excitations and derive a necessary and sufficient condition for quantum fractal liquids to be free from string-like logical 
operators. Several examples of quantum fractal liquids arc also studied. In particular, the Cubic code is shown to be 
unitarily equivalent to a certain model of quantum fractal liquids, which we call a second-order model. 

In section |Vj we compute topological entanglement entropies of quantum fractal liquids to establish connection 
between ground state and excitation properties. We demonstrate that topological entanglement entropies defined for 
tube-like regions can distinguish quantum fractal liquids from conventional topologically ordered systems. We give 
a physical interpretation based on propagations of quasi-particle excitations and their quantum dimension. Another 
important open question concerns coding properties of quantum fractal liquids. In section |VI[ we discuss coding 
properties of quantum fractal liquids by developing a mathematical method of finding lower bounds on code distances. 
We argue that quantum fractal liquids saturate the theoretical limit on quantum information capacity in three spatial 
dimensions, but without proof. Section [VII| is devoted to conclusions and open questions. Technical tools to analyze 
a family of translation symmetric stabilizer Haniiltonians are presented in appendix [X] which may be of independent 
interest. 

Some comments on the paper follow. We adopt the stability against local perturbations as the definition of 
topological order. By TQFT, we mean an axiomatic formulation by Atiyah which admits only a finite number of 
degenerate ground states [H]. For a review of TQFT, see [H]. By topological spin liquids, we mean gapped spin 
systems without local symmetries, i.e. topologically ordered spin systems. A connection between topological order 
and discrete gauge theory was pointed out by Read and Sachdev [7] and Wen [S] . Discussion on symmetry-protected 
topological order is beyond the scope of this paper. Discussion on gapless quantum spin liquids is beyond the scope of 
this paper. Our construction of quantum fractal liquids is theoretically motivated, and its relevance to experimental 
realization may not be immediately clear. Some technical tools are borrowed from a recent work by Haah |25j which 
aims to classify translation symmetric stabilizer Haniiltonians via commutative algebra. See [20! for a few other 
examples of string-free models proposed by Kim. See |21j for a model with pairs of fractal and string-like logical 
operators proposed by Castelnovo and Chamon. Limit cycles in RG transformations are not new phenomena as 
observed in the Efimov effect in non-relativistic conformal field theory with imaginary scaling dimensions [261 - 128] . As 
for relativistic field theory, it remains open for Z) > 2 [521 130] • See [3T] and subsequent papers by others for a recent 
proposal and debate. Results on coding properties of quantum fractal liquids, in section [VTl are based on an unproven 
conjecture. 
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II. TOPOLOGICAL SPIN LIQUID 



A. Topological spin liquid and string-nets 



In this section, we review how condensation of extended objects emerge in topological phases and present a physical 
picture of quantum fractal liquids, condensation of fractal objects. We begin with the simplest string-net model, 
known as Z2 spin liquid (or the Toric code) [HEJ] (Fig[T|a)). Consider a square lattice where qubits live on edges of 
the lattice with periodic boundary conditions. The Hamiltonian is 

s p rGs rGp 

where s represents a star and p represents a plaquette. Pauli X and Z operators act on each qubit as follows: 
Z\0) — |0), Z\l) = — 11), X\0) — |1) and X\l) = |0). The model is exactly solvable as interaction terms As and Bp 
commute with each other, and ground states satisfy 

As\^) = IV), Bp\^) = \^), ys,p. 

A ground state can be viewed as condensation of string-like extended objects. Consider a trivial product state jO)*^^ 
over the entire lattice {N is the total number of qubits) and observe that Bp\Q)'^^ = jO)**^. The following is a ground 
state: 



(1) 



since As{l + Ag) = 1 + Ag. The normalization factor is omitted. It is a superposition of terms As-^As^As^ ■ ■ ■ |0)®^. 
Since As is a product of Pauli-X operators, it flips the signs of qubits: |0) |1). Then a term As\0)'^^ can be viewed 
as a state with one small loop on a "dual lattice" , and a term Ag-^Ag^ 10)**^ with neighboring stars si and S2 is a state 
with a larger loop (Fig. [ijb)). In general, Ag-^Ag^Ag^ ■ ■ ■ jO)**^ is a state with loops of various sizes and shapes. A 
ground state is a superposition of all the loop states (Fig. [TJc)): 



'loop, 



V7 



(2) 



where 7 represents an arbitrary loop configuration. Therefore a ground state is condensation of fluctuating string-like 
objects with Z2 gauge symmetry. Note the model has four-fold degeneracy on a torus. 
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FIG. 1: Z2 spin liquid (the Toric code), (a) The Hamiltonian. (b) Loop states on a dual lattice, (c) Condensation of loops, 
(d) Logical operators. 
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One can construct a general quantum many-body system with several types of strings constrained by discrete gauge 
symmetry. Levin and Wen derived the most general form of wave-functions that are represented as condensation 
of string-like extended objects on a two-dimensional lattice by further assuming that wave-functions possess scale 
invariance and correspond to fixed-points of RG transformations |16j . Indeed, a ground state of Z2 spin liquid has 
scale invariance as it is a superposition of loops of all the different sizes and shapes. Note that scale invariance is 
required for systems described by TQFT since they must be invariant under continuous deformations. Yet, scale 
invariance is not a necessary condition for the presence of topological order. As we will see, quantum fractal liquids do 
not have full continuous scale symmetries. Instead, they have discrete scale symmetries where systems are invariant 
only under a limited set of scale transformations and ground states correspond to limit cycles of RG transformations. 

Geometric properties of extended-objects can be characterized by topological properties of global symmetry opera- 
tors. Formally, symmetries of the Hamiltonian can be captured by unitary transformations that leave the Hamiltonian 
invariant: 

U^'HU = H => [H,U]= 0. (3) 

Interaction terms and Bp are gauge symmetry operators for Z2 spin liquid since [j4s,iJ] = [Bp^H] — where the 
ground state space is an invariant subspace under actions of interaction terms. There also exist non-trivial symmetry 
operators which act non-trivially inside the ground state space (see Figjljd)): 

which have non-trivial winding on a torus. 

Since non-trivial symmetry operators commute with the Hamiltonian, they do not change the energy of the system. 
Yet, they cannot be written as products of Ag or Bp and transform degenerate ground states into each other. Recall 
that lipioop) is a condensation of loops that can be shrunk into a vacuum under Z2 gauge symmetry. An application 

(X) (X) 

of £q to IV'ioop) creates a non-trivial loop winding around the lattice, and \4'ioop} is condensation of loops with 

(X) 

non-trivial winding in the x direction. Similarly, £\ lipioop) is condensation of loops with non-trivial winding in the 
y direction. Four degenerate ground states may be indexed by winding numbers as \jx) ® \ly) with ^xily = 0,1 

where and 7^ represent the presence or absence of windings in the x and y directions respectively. Then non-trivial 

( x^ ( x^ 

symmetry operators and (\ act like Pauli-X operators on a pair of logical qubits \jx) It^)- It is convenient 
to represent their commutation relation as follows 




where operators in the same column anti-commute with each other while logical operators in different columns 
commute with each other. Anti-commuting pairs of non-trivial symmetry operators can be viewed as logical Pauli- 
Z and Pauli-A" acting on logical qubits \^x) "8) lly)- These non-trivial symmetry operators are often called logical 
operators in quantum information community as they can rewrite encoded logical qubits in the ground state space. 
Note Z2 spin liquid is a good quantum error-correcting code as only global operators with non-trivial winding can 
change encoded logical qubits, and no local errors can destroy stored quantum information. This is the origin of the 
stability against local perturbations in topological phases from quantum information perspective [2 133) . 

Another useful way of revealing topological properties of Z2 spin liquid is to consider its one- dimensional limit. 
Let us shrink the y axis of the lattice and consider a one-dimensional limit. Then a system consists of a pair of 
non-interacting ferromagnets in the x direction. Similarly, one obtains another pair of non-interacting ferromagnets 
by shrinking the x axis. Therefore, two-dimensional Z2 spin liquid can be constructed by weaving one-dimensional 
ferromagnets coherently such that its one-dimensional holographic images are ferromagnets. 

One can generalize Z2 spin liquid to higher-dimensional systems. For instance, the D-dimensional Toric code has 
anti-commuting pairs of m-dimensional and {D — rn)-dimensional logical operators: 

m-dim — m)-dim m : integer 

Its ground states can be viewed as condensation of m-dimensional extended-objects, or (£) — TO)-dimensional extended 
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objects in a dual description. In general, for a quantum many-body system described by TQFT, extended objects 
(Wilson loops or surfaces) have continuous geometries with integer dimensionality. In fact, the dimensional duality 
of non-trivial symmetry operators is a consequence of the Poincare duality for systems described by TQFT and can 
be derived analytically from continuous deformability of logical operators as shown in |34j . 



B. Emergence of fractal geometry 



In conventional topologically ordered systems, extended objects have continuous geometries with integer dimen- 
sionality. Yet, one may construct a quantum many-body system with fractal objects which have discrete geometries 
with non-integer dimensionality. In this subsection, we give a physical picture of quantum fractal liquids by reviewing 
how fractal geometries arise in classical spin systems. 

Consider a square lattice where L x L spins live on vertices of the lattice. The Hamiltonian is 

L-2L-1 , ^ ^^^^ 

where we represented interaction terms graphically as a matrix. We denote spin values at j) for i, j = 0, • • • , L — 1 
as Sij = 0, 1. Ground states must satisfy Ilijip = tjj, and 



(mod 2) 



(4) 



for all i and for j = 0, • • • , L — 2. The following is a ground state: 



1, 0, 0, 0, 0, 0, 0, 0, 

1, 1, 0, 0, 0, 0, 0, 0, 

1, 0, 1, 0, 0, 0, 0, 0, 

1, 1, 1, 1, 0, 0, 0, 0, 

1, 0, 0, 0, 1, 0, 0, 0, 

1, 1, 0, 0, 1, 1, 0, 0, 

1, 0, 1, 0, 1, 0, 1, 0, 

1, 1, 1, 1, 1, 1, 1, 1, 



where the upper-left corder corresponds to so,o- Configuration of sites with spin value 1 forms the Sierpinski trian- 
gle whose fractal dimension is log 3/ log 2 (see Fig. |4]^a)). It is interesting to observe that a ground state does not 
have translation symmetries while the parent Hamiltonian is translation symmetric. We call this spontaneous break- 
ing of translation symmetries strong-breaking of translation symmetries in contrast to weak-breaking of translation 
symmetries introduced in 

The model has a large number of degenerate ground states. Let us pick up arbitrary spin values on the first row 
of the lattice s = (so,Oi si^Oj S2,0j • • • ,sl~i,q)- Then, according to Eq. Q, spin values on other rows of the lattice 
are determined. Since there are 2^ possible choices for s, there are 2^ degenerate ground states. It is convenient 
to view the model as a time-evolution of one-dimensional cellular automaton where spin values on lower rows are 
computed via an update rule in Eq. Q. It is well known that one-dimensional cellular automata with linear update 
rules generate a variety of fractal geometries (see for a review). One may consider a general class of classical 
spin models with fractal ground states by designing interaction terms which imitate update rules of one-dimensional 
cellular automata. Such classical fractal liquid is particularly useful for storing classical bits of information securely 
as it asymptotically saturate a theoretical limit of classical information storage capacity of gapped spin systems ['I'd] . 

The Sierpinski triangle model has liquid-like order, but is different from conventional classical spin liquids such 
as anti-ferromagnetic Ising models on geometrically frustrated lattices |36H40) . Due to unconventional three-body 
interactions, the model does not have magnetic order at any temperature including T — 0. A zero temperature 
thermodynamic entropy is large, but not extensive: S = VN. In the Sierpinski triangle model, whether two spins 
have ferromagnetic or anti-ferromagnetic interactions depends on a spin value of the third spin, which leads to 
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emergence of peculiarly regulated fractal geometries. Weather the model may select an ordered ground state via order 
by disorder mechanism is not known. In this paper, we refer to a family of the Sierpinski triangle model as classical 
fractal liquids despite technical subtleties mentioned above. 

The model does not have power law decay of two-point correlation functions as observed in conventional classical 
spin liquids. Instead, we point out that it has an oscillatory power-law behavior with imaginary scaling dimensions, 
exhibiting discrete scale symmetries. Consider the following three-point correlation function: 



(5) 



In the ground space manifold, C(r) is 



C{r) = 1, r 2™ 
= 0, r^ 2" 

with oscillatory behaviors in log?', instead of r. At finite temperature, the three-point correlation function reads 

log 3/ log 2 



C{t) = (1 - 2pY 



with p 



= ^ / {e ^ + e^) where the exponent has a peculiar dependence on a fractal dimension. The correlation 



function for all r is 



C{r) oc exp(- A • r'°s 3/ log 2) . ^ 



27rfc 

r iog2 



k— — c 



where oscillatory behaviors in log r are represented by power law with imaginary scaling dimensions, which is charac- 
teristic of systems with discrete scale symmetries [27] . 



C(r) 



1 2 
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FIG. 2: Discrete scale symmetries and imaginary scaling dimensions in three-point correlation function. 

Geometric properties of degenerate ground states in the Sierpinski triangle model can be captured by geometric 
shapes of logical operators: 



Z, I, I, I, I, I, I, I, 

I, I, I, I, I, I, I, I, 

I, I, I, I, I, I, I, I, 

I, /, /, /, /, /, /, /, 

/, /, /, /, /, /, /, /, 

/, /, /, /, /, /, /, /, 

/, /, /, /, /, /, /, /, 

/, /, /, /, /, /, /, /, 



X, /, /, /, /, /, /, /, 

X, X, /, /, /, /, /, /, 

X, I, X, /, /, /, /, /, 

X, X, X, X, I, I, I, I, 

X, I, I, I, X, I, I, /, 

X, X, I, I, X, X, I, I, 

X, I, X, I, X, I, X, I, 

x^ x^ x^ x^ x^ x^ x^ x^ 



where they commute with all the interaction terms, but cannot be represented as products of interaction terms. So, 
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it has the foUowing pairs of logical operators: 

0-dim O l^-dim. 

log 2 

While the model has a fractal logical operator, a ground state does not have any quantum fluctuation, and is not 
topologically ordered since its partner is a trivial logical operator with zero-dimensional geometry. To have topological 
order, both logical operators must have topologically non-trivial geometries (i.e. they must be finite-dimensional). 

Given the large scope of theory of fractal geometries, it is natural to wonder if fluctuations of fractal objects can 
describe general topologically ordered quantum systems that are beyond descriptions of TQFT. In the reminder of the 
paper, we present a large class of quantum spin systems with fractal liquid order and analyze its physical properties. In 
particular, we construct a model of quantum spin systems which has pairs of anti-commuting fractal logical operators: 

fractal-dim <-> fractal-dim 

by combining features of both Z2 spin liquid and classical fractal liquids. 

C. Topological phase transition 

A quantum many-body system with topological order can be viewed as condensation of extended objects with a 
variety of geometric shapes. It is natural to expect that two ground states with different types of extended objects 
exhibit totally different physical properties, and thus belong to distinct topological phases. In this subsection, we make 
this intuition more precise by arguing that two spin systems with topologically different classes of logical operators 
are always separated by quantum phase transitions. 

We begin by reviewing a basic idea of classification of quantum phases by following [T71 ST] . The notion of quantum 
phases characterizes long-range ground state properties of a many-body quantum system with mass gap at zero 
temperature. Ground states in different quantum phases cannot be connected smoothly at the thermodynamic limit. 
Let us consider two ground states {ipA) and |'0b) of two different gapped Hamiltonian Ha and Hb and ask if they are 
separated by quantum phase transitions (non-analytic changes of ground state properties) or not. Two ground states 
IV'a) and IV's) are said to be in different quantum phases when there always exist some quantum phase transitions 
between and Hg regardless of choices of paths connecting and Hb- Conversely, if there exists a continuous 
change from Ha to Hb without crossing quantum phase transitions, two ground states {i/ja) and \iPb) are in the same 
quantum phase. In summary, one has the following classification principle: 

Being always separated by QPTs <^ Different quantum phases 

An equivalent, but more convenient way of classifying quantum phases uses local unitary transformations. Two 
ground states \iPa) and lips) are considered to be in the same quantum phase when there exists some local unitary 
transformation connecting \4'a) and \4'b)- By local unitary transformations, we mean transformations generated by 
a set of geometrically local quantum operations, applied for a finite duration. On the other hand, when there is no 
local unitary transformation connecting \tjjA) and \'(Pb), they are in different quantum phases: 

Being connected only by global unitary <^ Different quantum phases 

Note that only global unitary transformations can change long-range properties of ground states. 

These two classification principles of quantum phases, namely quantum phase transitions and local unitary transfor- 
mations, are equivalent under appropriate assumptions. If two gapped Hamiltonians Ha and Hb can be transformed 
into each other continuously without closing the energy gap, correlation lengths of ground states remain finite, and 
ground states at each stage of transformation can be approximated by applying some quasi-local unitary transforma- 
tions to original ground states [33]. Conversely, if \tpA) and \iPb) are connected by some local unitary transformations, 
one can always continuously transform Ha into Hb- Discussion so far is summarized in Fig. [3](a). 

The classification of quantum phases, based on continuous deformability of ground state wave-functions, reminds 
us of the study of topology in mathematics, which aims to classify geometric shapes of object based on continuous 
deformability. Roughly speaking, two objects are considered to be equivalent when they can be transformed into each 
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FIG. 3: Classification of quantum pliases. (a) Continuous deformability of ground states, (b) Continuous deformability of 
logical operators. 

other via continuous deformations (diffeomorphism) . Yet, if one cannot continuously deform an object to the other, 
they are considered to be topologically different. The similarity between classifications of quantum phases, based on 
continuous deformability of wave functions, and classifications of geometric shapes, based on continuous deformability 
of geometric objects, allows us to use the notion of topology in classifying quantum phases. Indeed, the following 
relation holds 

Logical operators are topologically different ^ Two systems belong to different quantum phases. 

The argument roughly goes as follows. Consider two systems with topologically distinct logical operators i and t' . 
Let us suppose that they belong to the same quantum phase. Then, there must be some local unitary transformation 
U such that UlW = i' . Yet, this is not possible since local unitary transformation can change geometric shapes of 
logical operators only continuously at the thermodynamic limit. Therefore, models with topologically different types 
of logical operators belong to different quantum phases and are always separated by quantum phase transitions. Note 
that there is no local unitary that transforms a string-like logical operator to a fractal logical operator, and thus, 
fractal models are different from conventional topologically ordered systems. "-4=:" of the above relation is proven only 
for stabilizer Hamiltonians with continuous scale symmetries (i.e. within TQFT) |43) . 

In summary, we expect that there will be four classes of quantum phases arising in stabilizer Hamiltonians. 

• Ferromagnetic phase has zero-dimensional logical operators and D-dimensional logical operators as in fer- 
romagnet. It possesses continuous scale symmetries with finite ground state degeneracy, and its ground states 
correspond to fixed points of RG transformations. Translation symmetries are not strongly broken. 

• Classical fractal phase has zero-dimensional logical operators and fractal logical operators as in classical 
fractal liquids. It possesses discrete scale symmetries with increasing ground state degeneracy, and its ground 
states correspond to limit cycles of RG transformations. Translation symmetries are strongly broken. 

• Topological phase has m-dimensional logical operators and {D — m)-dimensional logical operators where 
m > is an integer as in Z2 spin liquid. It is topologically ordered with stability against local perturbations. It 
possesses continuous scale symmetries with finite ground state degeneracy, and its ground states correspond to 
fixed-points of RG transformations. Translation symmetries are not strongly broken in the ground state space. 

• Quantum fractal phase has pairs of fractal logical operators as in quantum fractal liquids. It is topologically 
ordered with stability against local perturbations. It possesses discrete scale symmetries with increasing ground 
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state degeneracy, and its ground states correspond to limit cycles of RG transformations. Translation symmetries 
are strongly broken. 



III. CLASSICAL FRACTAL LIQUID 
A. Fractal and algebraic symmetry 



(a: 



(b) 



(c) 



FIG. 4: Examples of fractal geometries generated by polynomials. Unfilled elements represent zero, gray elements represent one, 
and red elements represent two. (a) The Sierpinski triangle from f = 1 + x over F2 with the fractal dimension D — log 3/ log 2. 
(b) The Fibonacci model from f = 1 + x + x'^ over F2 with the fractal dimension V — log(l + \/5)/log2. (c) The generalized 
Sierpinski triangle from f = 1 + x over F3 with the fractal dimension D — log 6/ log 3. 

In this section, we present a theoretical framework to construct a family of classical fractal liquids. It is well known 
that a variety of fractal geometries can be generated by polynomials over finite fields. We begin with polynomial 
representation of the Sierpinski triangle (Fig. |4]ja)). Consider a polynomial f = 1 + x over F2 and its j-th powers: 



f 


= 1 






= 1 - 


\-x 


f 


= 1 - 


Vx^ 


f 


= 1 - 


Vx- 




= 1 - 


Vx^ 


f 


= 1 - 


Vx- 



where coefficients are computed modulo 2. More graphically, one has 
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where "+" signs are omitted. One will notice that the Sierpinski triangle emerges in a geometric pattern of non-zero 
coefficients in P (j = 0, 1, • • • ). 

The entire Sierpinski triangle can be represented as a single polynomial with x and y: 



fix,y) = l + fy + py^ + fy^ + 
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where j-th row is indexed by . More graphically, one has 
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where non-zero coefficients of x''y^ correspond to filled elements of the Sierpinski triangle at Note that row 

vectors of the Sierpinski triangle evolves in the y direction via applications of generating polynomial /. 
Another interesting example of fractal geometries is generated by / = 1 + a; + a;^ over Fj,: 
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Again, the entire fractal geometry can be represented as f(x,y) ~ 1 + /y + py^ + py^ + • • • . The model is often 
called the Fibonacci model since its fractal dimension is given by '^°\lf^ (Fig- 4j[b)). The above construction of 
fractal geometries can be generalized to polynomials over Fp (p > 2) with an arbitrary prime p. For instance, consider 
J ^ 1 + X over F3 (p = 3): 
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which is a generalization of the Sierpinski triangle for three-dimensional spins (Fig.[4jc)). 

The self-similarity in fractal geometries arises from discrete scale symmetries of generating polynomials. Consider 
an arbitrary polynomial / over Fp with prime p: 



/ Co -(- CiX -I- C2X 



C3X 



(6) 



where Cj = 0, • • • ,p — 1. Then its p-th power is 

= Co + CIX'' + C2X2P + C3X3P + • • • . (7) 

For instance, with / = 1 + 2x + x^ over F3, one finds 

/ = 1 + 2x + x^ /3 = 1 + 2x3 + x*^ p = l + 2x^ + x^s /^^ = 1 + 2x2^ + x^'^. 

Due to discrete scale symmetries, the generated fractal geometry f(x, y) = l+fy+Py^+Py'^ + - ■ ■ has a self-similarity 
where the same pattern appear repeatedly at different length scales. 

Fractal geometries do not possess gauge symmetries since growth of filled elements violates charge conservation 
where a single element may evolve into multiple elements of the same type in the y direction. This is in strong 
contrast with the fact that continuous geometries often have physical interpretations based on conservation laws 
associated with underlying gauge symmetries as in the case of TQFT [16j . Charge conservation in scale invariant spin 
models originate from group theoretical constraints imposed on the parent Hamiltonian. Fractal geometries obey a 
more general form of symmetries, which will be referred to as algebraic symmetries, for a possible relation to theory 
of algebraic geometry which concerns geometric structures of solutions of polynomial equations. 
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B. Polynomial representation of Pauli operators 

To construct parent Hamiltonians of classical fractal liquids, it is convenient to represent interaction terms by 
polynomials too. We begin with construction of Pauli operators for a system of qubits on a one-dimensional infinite 
chain, following |25j . Consider a polynomial / over F2: 

CO 

f= J2 ^J^'' c, -0,l. (8) 

j=-oo 

We define the corresponding Pauli operators as follows 

00 00 

z{f)= n ^(/)= n (9) 
j— — 00 j— — 00 

where Zj and Xj are Pauli operators acting on j-th qubit. So, a polynomial / encodes positions of qubits where Pauli 
operators Zj or Xj may act. For instance, f = 1 + x + x'^ and Z{f) — Z0Z1Z2. 

The polynomial representation of Pauli operators is particularly useful for studying spin systems with translation 
symmetries since translations can be concisely described in terms of polynomials. For instance, consider a Pauli 
operator Z(f) = ZqZiZ2 for f — 1 + x + x^. Then, its translation in the i+ direction is given by Z1Z2ZJ,, whose 
polynomial representation is Z(xf): 

f = 1 + x + x'^ xf = x + x"^ + x^ 

Z{f) = ^ Zixf) = Z1Z2Z3. 

In general, Z{xf) is a translation of Z{f) in the i+ direction. Similarly, a translation in the i_ direction is given by 
Z{x~^ f). One may generalize this formalism to higher-dimensional systems by adding extra variables y, ■ ■ ■ . 
To gain more insights, let us represent a one-dimensional ferromagnet by polynomials over F2: 

H = -Y,Z{x^{l + x)) 

3 

where Z(x^{l + x)) — ZjZj+i. The Sierpinski triangle model is 

if = -^Z(a;V(l + 2: + a;z/)) 

ij 

where interaction terms are translations of Z{l+x+xy). In general, one may consider a classical translation symmetric 
Hamiltonian 

H = -Y^ Z{x'y^ • • • a) (10) 

hi,- 

with an arbitrary polynomial a{x,y, ■ ■ ■) where interaction terms are translations of Z(a). Ground states of a 
Hamiltonian obey 

z(xV.--«)l^) = IV') Vz,j,---. (11) 

The polynomial representation of Pauli operators becomes particularly powerful in analyzing commutation relations 
between Pauli operators .^(/) and X(^g). Since we are interested in translation symmetric systems, we would like to 
obtain commutation relations between .^(/) and translations of X((j). Let us imagine that we check commutation 
relations between Z(/) and X(x^ g) for all j and assign integers dj = 0, 1 as follows 

d,=0 for [Z{f),X{x^g)]^0 
= 1 for {Z{f),X{x^g)} = 0. 
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Based on dj, we define the commutation polynomial P{f,g) as follows: 

P{f,g)^Y.'^,x^ s.t Z{f)X{x^g)^{^lf'X{x^g)Z{f). (12) 







Thus, the commutation polynomial P{f,g) is a collection of commutation relations between Z{f) and X(x^g). For 
instance, with / = 1 + x + x'^ and g = 1 + a;, Z{f) anti-commutcs only with X{x^^g) and X(x^g). So, the commutation 
polynomial is P{f,g) = x~^ + x^. 

The commutation polynomial P{f,g) can be concisely written by introducing the notion of dual: 



CO 



^ c,x^ ^ ^ c,x-^ (13) 



J — — 00 J — — 00 

„-l 



where the dual / is obtained by taking x ^ x . Then, the commutation polynomial is given by the convolution 

P{f,9)^.f9 (14) 

For instance, one has fg — {1 + x + x'^){l + x~^) = x^^ + 2 + 2x + x"^ — x^^ + x^ for the above example. The proof 
of Eq. ( 14 1 is straightforward. Generalization to polynomials over Fp is straightforward. 



C. Classical fractal liquid 

We present general construction of classical fractal liquids. Consider a two-dimensional square lattice with L x i 
spins [L = 2"*). The Hamiltonian is 

H = -Y,Z{x'y'a), a = l + f{x)y (15) 

where f{x) is an arbitrary polynomial over F2 with x only. We put periodic boundary conditions both in the x and 
y directions, and assume that f{x) is reversible. (See appendix |A] for discussion on reversibility). 

In finding ground states of the Hamiltonian, it is convenient to find its logical operators. Z-type logical operators 
are trivial single Pauli operators ij = Zjq = Z{x^) for j = 0, • • • ■,L — 1 while X-type logical operators have fractal 
geometries: 

if^ = X{xH{x, y)), f(x, y) = l + fy+... + {fy)^-^ 

where f(x, y) is a fractal geometry generated by f{x). One can see that commute with all the stabilizer generators 
since the commutation polynomial between Z{1 + fy) and X{i{x, y)) is 

P(l + Jy,l + fy + ... + ifyf-') = (1 + fy){l + fy + ... + (/y)^-i) = 

when f{x) is reversible. We list all the logical operators as follows: 

f /(^) . . . /(^) 1 

l_ to ' '*L-1 ) 

where pairs of logical operators in the same column anti-commute with each other while pairs of logical operators in 
different columns commute with each other. So, there are k = L logical bits in total. 

(X) 

By using X-type logical operators £j , one can find all the ground states of a classical fractal liquid. Let us denote 
spin values at (i, j) as = 0, 1 for i, j = 0, ■ ■ ■ ,L — 1, and represent a ground state V' as 



= ^Sijx'^y^. 



(16) 
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Since the Hamiltonian consists only of Z-type Pauli operators, '0 = with = is a ground state of the Hamiltonian. 
(Recall Z\0) — |0) and Z\l) ~ — 11) in our notation). To find another ground state, one applies £q to = and 
obtains a fractal ground state: 

= l + fy + ... + {fy)^-^ = i{x, y) (17) 

One can find all the other ground states by applying fractal logical operators ^j"^"*- There are 2^ degenerate ground 
states, represented by 

0(7) = -fix)il + fy+... + ify)^-^) = 7(a:)f(x, y) (18) 

where ^{x) is an arbitrary polynomial with x only. Noting dim 7 — L, one finds k — L. 

Classical fractal liquids discussed so far are based on first-order cellular automata whose present states at t — t 
depend on states at < = r — 1. In higher-order cellular automata, the present states at t = r may depend on states at 
t = T — q, ■ ■ ■ ,T— 1 for q > I. One can construct classical fractal liquids based on higher-order cellular automata by 
taking 

a = 1 + hix)y + f2{x)y' + ■■■ + f,{x)y'^. 

However, it is generally difficult to write down spin configurations of higher-order classical fractal liquids explicitly. 

Since the model does not have gauge symmetries, its quasi-particle excitations violate charge conservation and 
propagate according to algebraic symmetries imposed by generating polynomial f{x). Recall that ground states 
satisfy Z{x'^y'-> a)il) = -0 for all and quasi-particle excitations may be viewed as violations of these algebraic 
constraints. It is convenient to represent positions of excitations by an excitation polynomial: 

E(x,y) = ^c.jx'y^ 

where an excited state ip' has 

Cij = Z{x'y^a)ip' = -H0' 
c,j = 1 Z{x'y^a)^' = -iP' 

such that a quasi-particle is present at if and only if Cij — 1. 

Excitations in classical spin liquids are caused by Pauli-X spin flips. Consider quasi-particle excitations caused by 
X{e{x, y)) where e(x, y) are polynomials representing positions of spin flips. Since anti-commutations between e(x, y) 
and Z{x^y^a) create quasi-particles at the excitation polynomial is 

E{x,y) = e{x,y)a. (19) 

For instance, if ^0,0 with e = 1 is applied, one has multiple excitations E{x, y) — a. Consider an isolated excitation at 
(0, 0). An application of Xq,o makes it propagate in the y direction to multiple excitations represented by f{x)y. So, 
quasi-particle excitations propagate via applications of f{x) just like a one-dimensional cellular automaton. Since the 
model does not have gauge symmetries, one cannot associate conserved charge to quasi-particle excitations. Indeed, a 
single quasi-particle may split into multiple quasi-particles of the same type. The energy barrier between degenerate 
ground states is 0{\og{L)) which leads to quasi-glassy thermodynamic relaxation as shown by Newman and Moore |22) . 
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' (D <D Z) fi^) 

FIG. 5: Propagation of quasi-particle excitations and violation of charge conservation. 



D. Limit cycle under RG transformation 



Classical fractal liquids have an interesting symmetry property concerning how their ground states behave under 
RG transformations. Let us consider the Sierpinski triangle model for L = 8. A ground state is 

1, 0, 0, 0, 0, 0, 0, o' 

1, 1, 0, 0, 0, 0, 0, 

1, 0, 1, 0, 0, 0, 0, 

1, 1, 1, 1, 0, 0, 0, 

1, 0, 0, 0, 1, 0, 0, 

1, 1, 0, 0, 1, 1, 0, 

1, 0, 1, 0, 1, 0, 1, 

1, 1, 1, 1, 1, 1, 1, 1 

RG transformations, denoted by RGy (i, j = 0, 1), pick up spins at (x, y) with x — i (mod 2) and y ~ j (mod 2): 



RGo.o(V') = 



1, 0, 0, 0' 
1, 1, 0, 
1, 0, 1, 

1, 1, 1, 1 

1, 0, 0, 0' 
1, 1, 0, 
1, 0, 1, 

1, 1, 1, 1 



RGo,i(V') = 



RGi.i(V') 



0, 0, 0, 

0, 0, 0, 

0, 0, 0, 

0, 0, 0, 

1, 0, 0, 0^ 
1, 1, 0, 
1, 0, 1, 

1, 1, 1, 1 



All the RG'ed states are ground states of the Hamiltonian for L = 4. Let us look at another ground state 

1, 1, 0, 0, 0, 0, 0, o' 

1, 0, 1, 0, 0, 0, 0, 

1, 1, 1, 1, 0, 0, 0, 

1, 0, 0, 0, 1, 0, 0, 

1, 1, 0, 0, 1, 1, 0, 

1, 0, 1, 0, 1, 0, 1, 

1, 1, 1, 1, 1, 1, 1, 1 

0, 0, 0, 0, 0, 0, 0, 
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and its RG transformations: 



RGo,o(V') 



RGi,o(V) = 



1, 0, 0, 0" 
1, 1, 0, 
1, 0, 1, 

1, 1, 1, 1 

1, 1, 0, 0" 
1, 0, 1, 

1, 1, 1, 1 

P, 0, 0, 



RGo,i(V') = 



RGi,i(^) = 



1, 0, 0, 0' 
1, 1, 0, 
1, 0, 1, 

1, 1, 1, 1 

0, 0, 0, 0' 

0, 0, 0, 

0, 0, 0, 

0, 0, 0, 



Again, RG'ed states are ground states of a smaller system. In general, one can define an arbitrary linear map from 
four spins in a 2 x 2 block to a single spin as an RG function. (Some non-linear mappings, such as a "majority vote", 
also work). One can see that RG transformations are stable against small noises added to ground state wave-functions. 
Ground states of the Sierpinski triangle model behave nicely under scale transformations by a factor of 2 only. If 
one performs a similar RG transformation by a factor of 3, RG'ed states are not ground states of the Hamiltonian 
anymore and flow to some disordered states. So, the model has scale symmetries under some limited set of scale 
transformations. This is a striking contrast with the fact that a ferromagnet, a spin model with continuous scale 
symmetries, look always the same under any scale transformations. Note that among four RG'ed ground states, two 
of them arc independent. 

It turns out that the presence of discrete scale symmetries is a general property of classical Hamiltonians with 
interaction terms Z{a) for an arbitrary polynomial a over Fp in any spatial dimensions. For simplicity of discussion, we 
keep concentrating on two-dimensional cases with a = l + f{x)y. Let us represent a ground state as V"™ = 'y{x)fm{x, y) 
where L = p"^ and f„i{x,y) is the polynomial representation of the fractal. Wc denote RG functions as KGij 
{i,j = 0, • • • ,p — 1) which pick up spins at {x,y) where x = i (mod p) and y = j (mod p). Then, RG'ed states 
RGij(^m) is always a ground state of the Hamiltonian for L = p"^~^. In particular, a polynomial 7' (a;) satisfying the 
following equation always exists: 



RGij{'y{x)im) = 7'(a;)fm_i. 



(20) 



The proof is immediate by recalling discrete scale symmetries of polynomials over Fp. For f = 1 + x over F3, one has 

1, 0, 0, 0, 0, 0, 0, 0, 0, 

1, 1, 0, 0, 0, 0, 0, 0, 0, 

1, 2, 1, 0, 0, 0, 0, 0, 0, 

1, 0, 0, 1, 0, 0, 0, 0, 0, 

1, 1, 0, 1, 1, 0, 0, 0, 0, 

Wl) = 1, 2, 1, 1, 2, 1, 0, 0, 0, 

1, 0, 0, 2, 0, 0, 1, 0, 0, 

1, 1, 0, 2, 2, 0, 1, 1, 0, 

1, 2, 1, 2, 1, 2, 1, 2, 1, 
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and 



RGo,o(^(l)) 



1, 0, 0, 0, 0, 0, 0, 0, 0, 

1, 1, 0, 0, 0, 0, 0, 0, 0, 

1, 2, 1, 0, 0, 0, 0, 0, 0, 

1, 0, 0, 1, 0, 0, 0, 0, 0, 

1, 1, 0, 1, 1, 0, 0, 0, 0, 

1, 2, 1, 1, 2, 1, 0, 0, 0, 

1, 0, 0, 2, 0, 0, 1, 0, 0, 

1, 1, 0, 2, 2, 0, 1, 1, 0, 

1, 2, 1, 2, 1, 2, 1, 2, 1, 



For the Fibonacci model / = 1 



over Ft , one has: 



^-(1) = 



1, 0, 0, 0, 0, 0, 0, 0, 0, 

1, 1, 1, 0, 0, 0, 0, 0, 0, 

1, 0, 1, 0, 1, 0, 0, 0, 0, 

1, 1, 0, 1, 0, 1, 1, 0, 0, 

1, 0, 0, 0, 1, 0, 0, 0, 1, 



2, 0, 0, 0, 0, 0, 0, 0, 0, 

2, 2, 0, 0, 0, 0, 0, 0, 0, 

2, 1, 2, 0, 0, 0, 0, 0, 0, 

2, 0, 0, 2, 0, 0, 0, 0, 0, 

2, 2, 0, 2, 2, 0, 0, 0, 0, 

2, 1, 2, 2, 1, 2, 0, 0, 0, 

2, 0, 0, 1, 0, 0, 2, 0, 0, 

2, 2, 0, 1, 1, 0, 2, 2, 0, 

2, 1, 2, 1, 2, 1, 2, 1, 2, 



and 



RGoo(^(l)) 



"1, 


0, 


0, 


0, 


0, 


0, 


0, 


0, 


0, 


1, 


1, 


1, 


0, 


0, 


0, 


0, 


0, 


0, 


1, 


0, 


1, 


0, 


1, 


0, 


0, 


0, 


0, 


1, 


1, 


0, 


1, 


0, 


1, 


1, 


0, 


0, 


1, 


0, 


0, 


0, 


1, 


0, 


0, 


0, 


1, 



RGoi(V(l)) 



1, 1, 0, 0, 0, 0, 0, 0, 0, 

1, 0, 0, 1, 0, 0, 0, 0, 0, 

1, 1, 1, 1, 1, 1, 0, 0, 0, 

1, 0, 1, 1, 1, 1, 0, 1, 0, 

1, 1, 0, 0, 1, 1, 0, 0, 1, 



Discrete scale symmetries provide an useful algorithm to compute the fractal dimension of f(x, y) We illustrate 
the algorithm for the Fibonacci model: / = 1 + x + over F2. We denote a ground state with an initial condition 7 
as ■0(7) • Then, renormalization of ground states i^mi^) and ■0m(H-x) gives the following ground states for L = 2™^^: 



RGoo(V^™(l)) =^^™-i(l) 
RGio(^™(l)) =i^™-i(0) 
RGoi(V^„,(l)) -7/',„_i(l + x) 
RGii(i^„,(l)) = i^,„_i(l) 



RGoo(^rn(l + a;)) =i/',„-i(l) 

RGio(V'm(l+-^)) =0m-l(l) 

RGoi(Vm(l + a;)) ^ i!„^-i{l) 

RGii(?/'m(l + X)) = llJm-l{x) 



Let us denote the weights of '!/'m(l) and '!/'m(l + x) as Aj^ and i?.,„. Then, one has 



B„ 



2 1 
4 



Am-l 
Brn-1 



(21) 



^ log I + 

This matrix has eigenvalues 1 ± v5, and thus. Am and Bm scale as 0{L ) for large L. 

The above RG transformations concern classical fractal liquids on a finite lattice. If one performs RG transforma- 
tions on an infinite lattice, RGij{tp) becomes a group operation where RGij{'ip) is a linear map inside the ground state 
space. In the case of a ferromagnet, the RG functions are always trivial; IlGij{ip) — ip since ^ is spatially uniform. 
Yet, for classical fractal liquids, translation symmetries of ground states are strongly broken, and as a result, RGij(V') 
may be different from ip in general. 

This gives an interesting possibility of limit cycle behaviors under RG transformations. Consider f = 1 + x over F3. 
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Let us apply an RG transformation for a ground state ■0(1) = ^{^tV) where i{x^y) = 1 + fy + py^ + • • • is defined 
on a infinite lattice. Then, RGi2('!/') gives the following sequence 

^(1)^^(2)^V(1)^V'(2)^--- 

where a ground state ?A(1) jumps to a different ground state 0(2), and the RG sequence exhibits a limit-cycle behavior. 
Next, for / = 1 + X + over F2, consider a ground state '0(1) — ^(2;, y)- Then, one has the following sequence under 
RGoi: 

V'(l) 0(1 + X) ^ 0(1) ^ 0^(1 + X) • • • 

which is also a limit cycle. 

Finally, consider f = 1 + x + over F5. We list some of its ground states as follows: 



1 

1 1 1 

1 2 3 2 1 

13 12 13 1 
.1 4 1 4 1 4 1 

1 3 

14 4 3 
10 4 12 3 

1 1 2 1 3 
1221233433 



1 1 

12 2 1 
1 3 3 1 
1 4 4 3 3 4 4 1 
1 04100140 1 

1 4 

10 4 
1114 4 4 
1 2 3 2 4 2 3 4 
1312434424 



1 2 

13 3 2 

1 4 2 3 2 

1 2 4 2 2 

1131142442 



A transformation RG02 generates the following limit cycle and fixed point: 

l-^l + a;^l + 3a;^l + 2x^l, l + 4a;^l + 4a; 

where ground states are represented by their initial conditions. A transformation RG04 leads to 

1 ^ 1 + .T ^ 1, 1 + 2a: ^ 1 + 4a; ^ 1 + 3.T ^ 1 + 3a;. 

These sequences are shown in Fig. |6] 




(b) 




P i+2x 



FIG. 6: Limit cycles in RG transformations for a classical fractal liquid with f = 1 + x + over F5. (a) RGo2. (b) RG04. 

How ground states transform under RG transformations can be formulated as a transition rule of polynomials. Let 
us define RG functions for one-dimensional polynomials, denoted by RGj as follows; for 7 = ^ . ax^, one has 

RG,(7) = E^f'+^-^'- (22) 

i 

In other words, we pick up coefficients of x^'^^^ . Let us consider a ground state ipij)- Then, its renormaHzation is 
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given by: 

RG,,(7/;(7))=^(7'), i' = RG,{jf). (23) 

We do not know whether hmit cycle behaviors lead to some interesting physical phenomena that are observable. 
Perhaps quantum criticality of classical fractal liquids (if exists) may have some reminiscent of limit cycles. We have 
considered wave-function renormalization so far. If one performs Hamiltonian renormalization, one obtains p copies 
of the original Hamiltonian after each transformation. 



IV. QUANTUM FRACTAL LIQUID 

A. Z2 spin liquid and polynomial 

In this section, we present a general framework to construct a family of quantum fractal liquids which are conden- 
sation of fractal objects. We begin by representing Z2 spin liquid by polynomials. Following |43j . we group two qubits 
into a single composite particle (Fig. [t]) such that composite particles live on vertices of a square lattice and obtain: 



qiZ) _ 



ZaZb, Za 
Zb, I 



I, Xa 
Xb, XaXb 



where each qubit inside a composite particle is labelled by A and B. In polynomial representation, the parent 
Hamiltonian is 



H 



^^UV(l + 2/)/ ^ 



where the upper (lower) row represents Pauli operators acting on A (B). Interaction terms arc translations of 

1 + 2: \ _ / H- 



1 + y 



X 



1 + a;-i 



In this form, it is immediate to see that Z2 spin liquid consists of one-dimensional ferromagnets at its one-dimensional 
limit since classical Hamiltonians generated by Z{1 + x) and Z(l +y) correspond to one-dimensional ferromagnets in 
the X and y directions respectively. 




X z 

o o 
xo ox ZO o 



I 



ZaZb 


Za 




Xa 


Zb 




Xb 


XaXb 



FIG. 7: Reduction of Z2 spin liquid. Two qubits are grouped into a composite particle. 
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Logical operators are 



Z 





i + x + x^ ^ — 





6P 



X 



It is worth representing them graphically as follows: 



'^0 



I, 



and 



(X) 



Xb, I, I, 

Xb, I, I, 

Xb, /, /, 

Xb, I, I, 



Zb, Zb 



I, I 



l + y + y^ + -- 






Za, I, I, 

Za, I, I, 

Za, I, I, 

Za, I, I, 



Xa, Xa, 
I, I, 



Xa, Xa 
I, I 



One can see that logical operators commute with interaction terms by computing commutation polynomials. For 

instance, = Z{0, 1 + a: + a;^ H )^ commutes with X{1 + 1 + x~^) since {1 + x + x"^ -\ )(1 + .x) = 0. 

One can generalize construction of Z2 spin liquid with arbitrary polynomials a{x, y,- ■ ■) and P{x, ?/,•••): 



Z 



X 



(24) 



Note that interaction terms commute 



where a and /3 are duals of a and /3 obtained by taking x ^ x ^ , y ^ y ^, 
with each other as their commutation polynomial is a/3 + /3a = over F2. A parent Hamiltonian is 



h=-j:z 



x^y^ 



■ a 



x^y^ ■ ■ ■ P 



x'^y^ ■ ■ ■ /3 
x^y^ ■ ■ ■ a 



Various physical properties of above Hamiltonians can be discussed from algebraic arguments on a and f3. See 
appendix [A| for further discussion. 



B. Quantum fractal liquid 

Consider a three-dimensional L x L x L square lattice where two qubits live on each site with L = 2™ and periodic 
boundary conditions. Quantum fractal liquids have 

a = l + fix)y, P = l+gix)z (25) 

where f{x) and g{x) are reversible polynomials over F2. More explicitly, interaction terms are translations of 

^f!+KV (26) 



l+g{x)zj' \l + f{x)y 

Therefore, interaction terms are characterized by a pair of fractal generators f{x) and g{x). In this sense, quantum 
fractal liquids can be viewed as a coherent combination of a pair of classical fractal liquids living on (a;,y)-plane and 
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(i, z)-plane respectively. Generalization to polynomials over ¥p is straightforward. 

Logical operators of quantum fractal liquids have fractal shapes which are generated by polynomials f{x) and g{x): 

fix,y) = l + fy + fy^ + ■■■ \x,y)^l + fy + Pf + ■ ■ ■ 
g{x, z) = l+gz + g^z^ H g{x, y) ^ 1 + gz + g^z^ H . 

Note i{x,y) lies on a (a;,y)-plane while g(x, z) lies on a (a;, z)-plane. Quantmir fractal liquids have k = 2L, and there 
are 2L of Z-type logical operators and 2L of X-type logical operators: 

-^[x%x,y))^ ) 

I j' -""[x^gix^z] J 

where i = 0, • • • , L — 1. Therefore, Z-type logical operators have geometric shapes of f(x, y) and g(a;, y) while X-type 
logical operators have geometric shapes of f{x,y) and g{x,y). 

To show that above operators are logical operators, we need to verify the following two things; a) they commute 
with interaction terms, b) they can be grouped into pairs of anti-commuting logical operators. One may see that 
logical operators commute with all the interaction terms by computing commutation polynomials. For instance, a 
commutation polynomial between i'jf '^ — Z{0, f{x, y))^ and stabilizers X(l+gy, 1+fx)'^ is given by {l+fy)t(x^ y) = 0. 
Logical operators obey the following commutation relations: 

f /(^) ... ... 1 

J '^O ' ' ^L-l' 'O ' ' I 

1 r^^) ... /-^) ... /(^) f ■ 

'O ' ' ^0 ' ' ) 

To see that quantum fractal liquids are topologically ordered, we begin by showing that they are good quantum 
error-correcting code with d — > oo for L ^ oo where d is the quantum code distance of the ground state space. A 
standard way to prove this considers a bi-partition of the system into two complementary subsets A and B and uses 
the following bi-partition formula which holds for arbitrary stabilizer codes |45j : 

gA + gs = 2k (27) 

where gA and gs represent the number of independent logical operator supported inside A and B respectively. Let 
us assume A to be a connected region with finite support. Then, its complementary subset B accommodates some 
(i, y)-plane and {x, £)-plane where all the 2k independent logical operators can be supported. So, one has gs — 2k. 
This leads to gA — 0. Therefore, weights of logical operators are not finite (i.e. unbounded), and d ~> oo for L — > oo. 
Another way uses a criteria presented in appendix [A| (see (201 also). Quantum fractal liquids are characterized by 
a pair of polynomials a = 1 + f{x)y and P = 1 + g{x)z which do not have common factors. Therefore, the system 
always has d — > oo for i — > oo (as long as /c ^ 0). 

For stabilizer Hamiltonians, being a quantum code (c? — >■ oo for L — >■ oo) automatically implies the presence of 
topological order with stability against local perturbations. In [33], Bravyi, Hastings and Michalakis proved that 
frustration-free Hamiltonians with an ability of quantum error-correcting code have stability against local perturba- 
tions when Hamiltonians satisfy a certain condition, called TQO-2. Roughly speaking, TQO-2 states that locally 
computed density matrices are consistent with ground states which are computed globally. One can check that quan- 
tum fractal liquids satisfy TQO-2, and thus have stability against local perturbations. Recall that quantum fractal 
liquids have 2^ ground states. Under a small but finite local perturbations, the energy splitting among these ground 
states is always exponentially suppressed, and the energy gap between the ground states and excited states remains 
finite. Note quantum codes with TQO-2 condition have efficient decoding algorithm based on RG-like transforma- 
tions ps] . 

We then discuss the number of degenerate ground states and its dependence on the system size. A key feature of 
quantum fractal liquids is that the number of logical qubits k has a fairly sensitive dependence on the system size 
L = {Li, L^). The number of logical qubits k is given by counting the dimension of solutions 7 satisfying the 
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following equation: 

/(a:)^^7(x) = g(x)^^7(x) = 7(x), = 1. (28) 

with k = 2dini7. For instance, with / = and g — x-> , one has k = 2 gcd{Li,iL2, j L^) where k depends crucially on 
the system size L — (Li, L2, -^s)- In general, it is a very challenging task to write down an explicit form of fc(Li, L2, L3) 
for a given pair of f{x) and g{x). Yet, k{Li^ L2, L^) has a nice symmetry property under scale transformations: 

k{pLi,pL2,pL3) ^ pk{Li,L2,Lz). (29) 

The proof is straightforward from discrete scale symmetries of polynomials over Fp. See appendix |A] for quantum 
fractal liquids with open boundary conditions. 

Ground states of quantum fractal liquids correspond to limit cycles under real-space RG transformations on an 
infinite lattice. To obtain RG transformations with meaningful attractors that do not flow to disordered states or 
trivial product states, one needs to apply some appropriate projection operators on sites that are to be coarse-grained. 
Consider a pair of qubits at site {i,j,£), denoted as \(f>)ije = \4'A)iji ® |0b)ij£, and apply the following projections to 
a ground state: 

{I + Zi® Z'g){I + ® X^b)\4>a).31 ® \<t>B)^ji- (30) 

Note that projection operators commute with each other, and projections are applied only to sites {i,j,£) with j ^ 
OT £ ^ modulo 2. As a result, pairs of qubits on sites {i, j, £) with j = £ = modulo 2 are completely decoupled from 
the rest where projectors are applied. With some speculation, one notices that stabilizer generators for remaining 
sites (z, j, £) with j = £ = modulo 2 are given by 

x(f:), .(3) 

and their translations that are generated by applications of x^* y^^ z^^ . (See for transformations of stabilizer 
generators under projections). This corresponds to two copies of original quantum fractal liquids. Let us pick up sites 
with (i, j', £) = (0, 0, 0) modulo 2 and throw away sites with (i, j, £) — (1,0,0) modulo 2 via some arbitrary projections. 
Taking x^ ^ x, ^ y and — )■ 2;, stabilizer generators are Z(a,/3)-^ and X{P^aY' . So, this RG transformation 
maps a ground state of a quantum fractal liquid into some ground state which may be different from the original. 
One can keep track of how ground states flow in RG transformations by looking at polynomial representation of 
fractal logical operators which can be analyzed in a way similar to ground states of classical fractal liquids. If one 
performs Hamiltonian renormalization, instead of wave-function renormalization, one obtains two copies of the original 
Hamiltonian after each transformation. For quantum fractal liquids with higher-order cellular automata rules, RG 
flows become more complicated. 



C. Quasi-particles and no-string rule 



We discuss properties of quasi-particle excitations in quantum fractal liquids. Without loss of generality, one 
can concentrate on excitations caused by Pauli-.Z type operators which flips X-type interaction terms. Following a 
treatment of classical fractal liquids, we represent positions of excitations as a polynomial: 



E{x,y,z) = 



lit 

CijiX y-'z 



where Cijg — 1 means an excitation is present at {i, j, £). We consider excitations caused by a Pauli operator Z{eA: es)^ 
which is given by a commutation polynomial between (e^, es)"^ and (/3, ay^ : 



E{x,y,z) = P 



1 + .92 



eA{l + gz)+eB{l + fy). 



(31) 
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Excitations caused by a Pauli operator Zg q q {ca = 1, es = 0) are given by 

E{x,y,z)^{l + gz) (32) 

while excitations caused by a Pauli operator Zq q q {ca = 0, = 1) are given by 

E{x,y,z) = {l + ,fy). (33) 

So, if an isolated excitation is present at (i, j, i) = (0, 0, 0), it will propagate to multiple excitations represented by gz 
via an application of Zq'^q (Fig- jsj^s-))- Similarly, it will propagate to multiple excitations represented by fy via an 

application of Zq^Jq (Fig. jsj^a)). In general, for a single isolated excitation, f{x) is applied when propagating in the 
y direction, and g{x) is applied when propagating in the z direction. Let us consider a one-dimensional excitation 
pattern, represented by e{x), located at j — £ = 0, and think of moving it to j = j' and i = (! while keeping the 
pattern of excitations one-dimensional. Then, the resulting excitation pattern is given by 

E[x,y,z) = e{x)!{xY g{xYy^' z'' . (34) 

This may viewed as time evolution of an initial condition e(a;), updated j' times by j(x) and times by g(x) 
respectively. 




FIG. 8: (a) Propagation of quasi-particles by /(a;) and gix) (b) A pair of localized excitations ei and 62 with elongated 
excitations e{ and 

Let us consider simple cases with /(x) — g[x) and consider propagations of an excitation originally at (0,0,0) to 
i = j' and £ — I' . Then, the resulting excitation pattern is given by 

E{x,y,z) = j{xY g{xYy^' z'' = ^{x^ y^' z'' . (35) 

So, if j' -|-^' = 0, excitations are single quasi-particles, and excitation energy remains finite. Therefore, quasi-particle 
excitations can propagate freely (without costing much energy) in the y — z direction which implies the presence of a 
string-like logical operator. Indeed, the following string-like operators are logical operators: 

Note that (y + z)^^^ is a string-like polynomial extending in the y — z direction. 

A naturally arising question is when quasi-particle excitations propagate freely and quantum fractal liquids have 
string-like logical operators. Without loss of generality, one can assume that f{x) and g{x) start from the origin, 
meaning that f(x) and g{x) have non-zero constant terms and have only positive powers. We say that f{x) and g{x) 
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are algebraically related when there exists some finite integers c/ and Cg such that 

fixp - g{xr^ (36) 

without considering periodic boundary conditions. Then, one notices that quasi-particles can propagate in the Cajj—CbZ 
direction. Therefore, there exist string-hke logical operators. 



It turns out that Eq. (36) is a necessary and sufficient condition for the presence of string-like logical operators: 



"No string- like logical operator" "Polynomials /(x) and 5(2;) are not algebraically related" . (37) 

To verify this equivalence relation, let us consider a pair of localized excitations ei{x,y,x) and e2{x,y,x) which are 
contained in a cubic box oi w x w x w sites with w L, which are separated by Ls- We assume Ls is sufficiently 
larger than w, and sufficiently smaller than L. We ask whether ei and 62 are connected by a string-like object or 
not. One can make quasi-particle excitations in ei{x,y,z) and e2{x,y,z) propagate by applying f{x) and g{x) such 
that they become some elongated line-like excitations whose lengths are 0{w) (see Fig.|8|b)). We denote polynomials 
corresponding to these elongated excitations as e[ and 63, and assume that e'^ is at j — £ — and 63 is at j = j' 
and 1 = 1' where |j'| -I- ~ OiL^). One may write e'j^ = e|(a;) and 62 = e*2[x)y^ . Then, one must have 
^2 = ^ifi^y 9{^Y in order to connect these excitations by some object. If a pair of excitations is connected by 
a string-like object, g^ must remain finite for large j' and £' . This requires / and g to be algebraically related; 
otherwise, one can prove that the size of grows at least linearly as \j'\ and grow. So, a pair of isolated 

excitations is not connected by a string-like object, and there is no string-like logical operators. One may construct 
a formal proof by following an argument in [57] . 

When there is no string-like logical operators, energy barrier separating different ground states is lower bounded 
by O(logL) as shown in [221 147] . The presence of large energy barrier is favorable for storing logical qubits securely 
in thermal environment. For sufficiently small system size L where entropic contributions are negligible, the memory 
time of encoded logical qubits can be estimated by an Arrhenius law r ^ exp(— A/T) where A is an energy barrier. 
For models without string-like logical operators, one has r ~ POLY(L) which is an improvement over r ~ 0(1) 
of conventional models with string-like logical operators. For large L, however, one does not see any significant 
improvement since entropic contributions dominate over energetic constraints as the number of ground states and 
local minima increases. With this subtlety, fractal models without string-like logical operators are said to be marginally 
self-correcting \W\ ■ This behavior is closely related to the thermal instability of fractal models as they undergo thermal 
phase transition at Tc = 0. So, the transition to quantum fractal glassy phase is a continuous one (i.e. a crossover), 
and is not sharp. See [48 [ 149] for proposals to overcome this challenge in quantum information community. 



D. Several examples 

Finally, we study several examples of quantum fractal liquids, 
(a) We begin with a trivial case: 



1 + y 
l + z 



over F2 



This is a stack of slices of two-dimensional Z2 model where each copy lives on a (y, z)-plane. It has pairs of string-like 
logical operators since / = 1 and .g = 1 are generators of strings. Similarly, for f — I and g = I over Fp {p > 2), the 
model is a stack of Zp spin liquids with p-dimensional qudits. 
(b) For f = X and 5 = 1, stabilizer generators are given by 

z(\+''y] over F2 



l + z 

This is a stack of slices of two-dimensional Z2 model, but each copy lives on a {x + y, z)-plane. It has string-like logical 
operators, running in the x + y direction. One can reduce this model to the model in (a) by a modular transformation 
X ^ X, xy ^ y and z — )■ z which corresponds to a lattice distortion. In general, two models are considered to be 
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equivalent when they are connected by modular transformations. 

(c) Let us look at an example with pairs of fractal logical operators and string-like logical operators: 

l + {l + x + x^)y\ 

, , over Fa, 

1 + z J 

the model has pairs of fractal logical operators and string-like logical operators. In particular, geometric shapes of 
fractal logical operators are generated by / = 1 + x + (the Fibonacci model) . Note that fractal logical operators 
live on a two-dimensional plane and string-like logical operators are perpendicular to the two-dimensional plane. 

(d) Some models do not have any logical qubits under periodic boundary conditions: 

Z \ , \ over F2 

\ l + z ) 

since / = 1+x is not reversible over F2. When the model is defined with open boundary conditions (as in appendix[A|, 
it has pairs of fractal logical operators and string-like logical operators. 

(e) Consider 

Z[ , ^ T ^ , over F2. 

\ 1 + (1 + x)z J 

One might think that the model has pairs of fractal logical operators, and is free from string-like logical operators. 
Yet, the model has hidden string-like logical operators running in the y — z direction. In fact, this model is unitarily 
equivalent to the following model 



l + {l + x)y 
1 + y-^z 



over F2. 



See appendix [A| for further discussion on local unitary transformations. 

(f) Let us consider the cases without any string- like logical operators for F2: 



l + {l + x + x'^)y 
1 + (1 + .T + x^)z 



Z[ , . ^ . JC. I over F2 



The model is free from string-like logical operators as / = 1 + a; + and g ^ 1 + x + x'^ are algebraically unrelated. 
Interaction terms are seven-body. 

(g) Next, let us consider a model over Fp (p > 2): 



1 + (1 + x)y 
1 + (1 + x'^)z 



Z\ -, , A ' 2;., ) over F3 



The model is also free from string-like logical operators as / = 1 + a; and g = 1 + are algebraically unrelated. 
Interaction terms are five-body. 

(h) The Cubic code can be viewed as a second-order quantum fractal liquid. In polynomial representation, one has 



\ + x -\- y + z 
1 + xy + yz + zx 



over F2 . 



The model can be mapped to the following second-order quantum fractal liquid through local unitary transformations 
and modular transformations: 



l + fix)y 
I + gi{x)z + g2ix)z'^ 



Z{ , , \ _ 2 I over F2 
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where 

f{x) = l + x + x'^, gi{x) = l + x, g2{x) ^ 1 + X + x'^ . 
See appendix [K\ for further discussion. 

V. TOPOLOGICAL ENTANGLEMENT ENTROPY 

In this section, we compute entanglement entropies of quantum fractal liquids in order to establish a connection 
between ground state and excitation properties. We demonstrate that topological entanglement entropies, defined for 
tube-like regions with holes, can distinguish quantum fractal liquids from conventional topologically ordered systems. 

We begin with a review of computational technique of entanglement entropies. For stabilizer Hamiltonians with 
topological order, calculation of entanglement entropies can be reduced to a simple counting argument. Consider an 
arbitrary stabilizer Hamiltonian, H — — Sj and a bi-partition into complementary subsets A and B = A. We 
denote a group of Pauli operators generated by Sj a.s S — {{Sj : Vj}) and its restriction into A as Sa- 

SA^iiUeS: U\b - /}) (38) 

where U\b is a restriction of U on B. Note that Sa consists of all the stabilizer operators supported inside A. If the 
stabilizer Hamiltonian possesses topological order, there is no logical operator supported inside A, and one has 

EA^\A\-dimSA (39) 

where Ea is an entanglement entropy, \A\ is the total number of qubits inside A, and dim5yi is the number of 
independent generators for Sa- 

Let us now compute an entanglement entropy of Z2 spin liquid, represented by Z[l + x,l + y)'^ , for a region A with 
wi X sites. One has \A\ = 2wiW2 and dim5 = 2(wi — l)(u'2 — 1), so Ea — 2(u;i + ^2) — 2. There is a constant 
correction term to the area-law of entanglement entropies. Yet, this correction is not universal since it is dependent 
on shapes of boundaries. To obtain a universal correction term correctly, one needs to cancel out non-universal 
contributions by adding or subtracting entanglement entropies for several regions. 



(a) 



A 


B 


C 




Bi 



Ci 



C2 



FIG. 9: Computations of topological entanglement entropies. 



There are two standard methods for this [50|, I51j . The first method considers the mutual information among three 
neighboring regions (Fig. [9]ja)): 

I{A, B, C)=Ea + Eb+Ec- Eab - Ebc - Eca + Eabc- (40) 

Since Ej^ = \R\ — dimiS^ for regions R, contribution from volumes \R\ exactly cancels. As for stabilizer generators, 
U e Sa implies U e Sab,Sca,Sabc- So, contribution from U G Sa cancels. Only stabilizers U e Sabc, but with 
U ^ SatSb,Sc,SabtSbctSca contribute to I{A^B^C). Therefore, one needs to count the number of independent 
stabilizer generators with supports on all A^ B and C. One may see that there is only one such stabilizer, so 
I{A, B, C) = —1. However, we note that there exist models of stabilizer Hamiltonians which are connected to product 
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states, but possess non-zero I{A, B,C). For instance, a model with a = 1 and /? = 1 + x + y + xy has I{A,B,C) = —2 
though it is connected to a product state through local unitary transformations. This discrepancy is perhaps because 
boundaries of topologically ordered systems are not necessarily described by conformal field theory on a precise level. 
The other method considers the following combinations of entanglement entropies (Fig. |9][b)) 



E* = Ea — Ebi — Eb2 + EciuC2 



(41) 



where A is a loop-like region. One needs to find stabilizer generators with U £ Sa and U ^ Sbi,Sb2,SciuC2- To find 
such stabilizers, it is convenient to find identity generators R (see appendix [A] for detailed discussion) : 

R^^x'y^, Ra = Rf3 = 



which creates an identity operator as Z{Ra, RfH)^ — I. Consider a subpart of R, given by R* — J27j=o^^y''- This 
creates the following loop-like stabilizer operator 



R*a 

R*P 



ZaZb 


Zb ■ 


■ Zb 


Za 


Za 


I ■ 


■ I 


Za 


Za 


I ■ 


■ I 


Za 


. Zb 


Zb ■ 


■ Zb 


I _ 



(42) 



which has supports on all the four regions _Bi. i?2, Ci, C2. One can construct a similar operator from X-type stabilizers. 
So, one has E* = 2. These stabilizer operators may be viewed as string-like extended objects. Since endpoints of loops 
correspond to quasi-particle excitations, these loop- like stabilizer operators characterize propagations of quasi-particles 
along the loop A. There are two independent quasi-particles circling around the loop in Z2 spin liquid. 

Let us proceed to discussion on topological entanglement entropies of quantum fractal liquids in three-dimensional 
systems. There are four possible types of topological phases; (a) the three-dimensional Toric code, (b) a stack of 
slices of two-dimensional Toric code, (c) quantum fractal liquids with string-like logical operators, and (d) quantum 
fractal liquids without string-like logical operators. The goal is to distinguish these four phases. For this purpose, we 
consider topological entanglement entropies defined for a tube-like region with a hole inside it (Fig. [Io| ) . We assume 
that the width of the tube w is sufficiently large, but is small compared to L. We add and subtract entanglement 
entropies in a way similar to Fig. |9]^b). We argue that, one can, in principle, distinguish four phases by measuring 
topological entanglement entropies for tube-like regions in various directions. 

(a) The three-dimensional Toric code is 





with three qubits per site. (See appendix [A| for higher-dimensional Toric code). One can form a loop-like Z-type 
stabilizer by multiplying Z{0, 1 + z, 1 + y)-^writh R* = J2j7=o V'' ^ tube-like X-type stabilizer by multiplying 

X(l -f x, 1 -f y, 1 -I- z)^ with R* — Yl,^=o Sjfco ^^U^z^- So, one has E* = 2 since there are only two independent 
quasi-particles. One obtains the same result for tube-like regions extending in the y oi z directions. So, topological 
entanglement entropies remain finite for tubes in any directions. This is related to the fact that propagations of 
quasi-particles are rotationally invariant at the thermodynamic limit and translation symmetries are not broken in 
any directions. 

(b) We then compute topological entanglement entropies for a stack of slices of two-dimensional Toric code. To be 
specific, we assume that each copy of the Toric code lies on a (J), z)-plane. For a tube-like region extending in the x 
direction, one has E* = 2L which corresponds to the number of quasi-particles that can circle around the tube. Yet, 
for a tube-like region extending in the y direction, one has E* = since quasi-particles are constrained to (y, £)-planes 
and they can propagate "along" the tube, without circling around the tube. This is related to the fact that translation 
symmetries are strongly broken in the x direction. 
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(c) Next, we compute topological entanglement entropies for quantum fractal liquids with strings. To be specific, 
we consider the case with f{x) ^ and g{x) = 1 where string- like logical operators are along the z direction. Identity 
generators are 

R = j{x)f{x,y){l + z + z^ + ■■■) 

where j{x) is an arbitrary polynomial with x. For a tube-like region extending in the x direction, one can form tube- 
like Z-type stabilizers with R* = 7(2;) J2i=o ^JJ=o ^^{fyY^^- There are L independent tube-like Z-type stabilizers, 
so one has E* = 2L. This corresponds to the total number of independent quasi-particles that can circle around the 
tube. In fact, one obtains the same result for a tube-like region extending in the y directions: E* = 2L. Yet, for a 
tube extending in the z direction, with some speculations, one has £■* = since quasi-particles can propagate only 
along the z direction, without circling around the tube. Note that translation symmetries are strongly broken in the 
X and y directions. 

(d) Finally, we discuss quantum fractal liquids without strings. Identity generators are 

R = 'y{x)f{x,y)g{x,z). 

For a tube-like region extending in the x direction, one can form L independent tube-like Z-types stabilizers, so 
E* = 2L. In fact, one obtains the same result for tubes in the y and z directions: E* = 2L. This is due to an 

observation that quasi-particles cannot propagate along the tube, regardless of its direction, and need to circle around 
the tube. Note that translation symmetries are broken in all the directions, yet, rotational invariance of topological 
entanglement entropies are persevered. 
Below, we summarize the results: 





Tube in x 


Tube in y 


Tube in 


3-dim Toric code 


0(1) 


0(1) 


0(1) 


2-dim Toric code 


0(L) 








Quantum fractal liquid (with string) 


0(L) 


0(L) 





Quantum fractal liquid (no string) 


0(L) 


0(L) 


0(L) 



One may notice that topological entanglement entropies become zero or finite when string-like logical operators are 

present along the axis of the tube. So, by measuring topological entanglement entropies for tube-like regions in various 
directions, one may distinguish quantum fractal liquids from conventional topologically ordered systems in principle. 
(There are some artificial counterexamples to our claim. For instance, three copies of stacks of two-dimensional Toric 
codes, put on three different directions, have 0{L) topological entanglement entropies in all the three directions. 
Stability of these fractal order parameter is an important problem, but is beyond the scope of this paper). Another 
useful method of detecting quantum fractal order is to see if translation symmetries are strongly broken in the ground 
state space or not. One may be able to deduce fractal generators f{x) and g{x) by analyzing behaviors of k{Li,L2, L^). 




FIG. 10: A tube-like region extending in the x direction. 



29 



VI. TOWARD SATURATION OF QUANTUM CODING CAPACITY 
A. Coding properties of classical fractal liquid 

In this section, we discuss coding properties of quantum fractal liquids and argue that they may saturate the 
theoretical limit on quantum information storage capacity. We begin by reviewing coding properties of classical 
fractal liquids [23]. Consider a two-dimensional classical spin system supported by local commuting frustration-free 
Hamiltonians: 

where Ilij- > and ground states satisfy = for all So, a ground state can be obtained by minimizing 

interaction terms independently and there is a finite mass gap between ground states and excited states. Such 
Hamiltonians are often called local codes as they can be used for reliably storing classical bits of information inside 
the ground state space. In particular, we think of encoding classical bits of information into degenerate ground states 
of a local code. The number of encoded logical bits and the classical code distance are denoted as k and d. Note that 
the code distance d is the minimal Hamming distance among 2^ degenerate ground states, and corresponds to the 
minimal number of bit flips that are needed to transform one ground state to another. So, the code distance d is a 
quantitative measure of reliability of encoding. One may naturally expect that it is difficult to encode a large amount 
of logical bits while keeping a code distance d large. This intuition on a tradeoff between k and d can be made more 
rigorous as summarized in the local code hound [52] : 

k^/d < 0{n) 

for two spatial dimensions. Note a proportional factor on the right-hand side depends on the range of interaction 
terms. 

Then, a naturally arising question is; "what kinds of local codes saturate the bound?" In |23j . we showed that some 
classical fractal liquids asymptotically saturate the local code bound: 

fc~0(i), dr^O{L^-') 

where e is an arbitrary small positive number. Note a code distance d is upper bounded by n by its definition, so such 
local codes are essentially the best local codes. 

Consider a two-dimensional classical fractal liquid generated by a polynomial / over Fp on an L x L lattice. We 
denote the classical code distance as c?„j for L — p'". From discrete scale symmetries of classical fractal liquids, one 
can rigorously prove the following: 

dm>L^, di=p^ (43) 

where I? < 2 is a fractal dimension of the code distance. Let us consider a generalized Sicrpinski triangle model with 
/ = |(a; -|- 1) for Fp. (A factor ^ is added to ensure /p = 1 for p > 2). Fractal dimension of a generalized Sierpinski 
triangle is 

log^^ 

— > 2 for p — > oo 

logp 

which monotonically increases and approaches to two. This implies that a Hamming distance between a trivial ground 
state and a fractal ground state f(a;, y) is 0{L?~'^) for large p. In fact, one can prove that the fractal dimension T) 
of the code distance also behaves as 



I? — > 2 for p ^ oo. 



(44) 
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Therefore, there exists a classical fractal liquid with 

where k^^^ is the number of encoded p-dimensional spins. So, it asymptotically saturated the local code bound. One 
may simulate p-dimensional spins with two-dimensional spins. Then one has 

nr^ 0{L^ -logp), fc-O(L-logp), dr^O{L^'') 

where the range of interaction terms is w ^ O(logp). 

A similar argument holds for higher-dimensional local codes too. In D spatial dimensions, the local code bound 
reads 

kdT> < 0{n). 

There exists a classical fractal liquid with 

which asymptotically saturates the local code bound. It is interesting to observe that coding properties of classical 
fractal liquids obey an area law where the number of encoded bits k behaves area-like while the code distance d 
behaves volume-like. Note that an original argument in [53] was applied to classical fractal liquids with open boundary 
conditions. By the use of discrete scale symmetries and some algebraic argument, one can analytically prove a similar 
result for classical fractal liquids with periodic boundary conditions. 

B. Coding properties of quantum fractal liquid 

We then discuss coding properties of quantum fractal liquids. For local quantum codes, supported by commuting 
frustration-free Hamiltonians, the following bounds are known to hold 52j : 

d<0{L^-^) (45) 
kdT^ < 0{n) (46) 

in D spatial dimensions where n = and d is a quantum code distance. If these bounds are tight, the best local 
codes would be 

k^O{L"-^), d-O(L^-i) 

where k is area-like and d is volume-like in D space-time dimension. For D — 2, the Toric code is known to saturate 
the bounds with k ^ 0(1) and d ^ 0{L). Yet, the problem of finding capacity saturating local quantum codes remains 
open for D >i. In particular, the previous best code (_D-dimensional generalization of the Toric code) has k ^ 0(1) 
and d ^ 0(LL^/^J) where [-D/2J is the largest integer which is smaller than (or equal to) D/2. A similar question for 
quantum low density parity check codes (quantum LDPC codes) remains open too 53J. Quantum LDPC codes are 
local codes constructed on a graph with low density of interaction terms, but without restriction of geometric locality 
of interaction terms. The previous best quantum LDPC code has d ~ 0{y/n\ogn) [Mj . 

Here, we argue that a family of three-dimensional quantum fractal liquids saturate the above bounds: 

fc-O(L), d^O{L^) 

based on an unproven conjecture. So, quantum fractal liquids may be the first examples which do better than the 
Toric code in higher-dimensional systems. 

If / and g are algebraically related, a system has string-like logical operators which leads to a small code distance. 
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The simplest choice of string-free quantum fractal liquids is 

+ g^-[x-2) (47) 

over Fp. Note that fP = gP = 1 under periodic boundary conditions. Here, we consider a model on a lattice with 
p X p X p sites. As proven in appendix]^ a lower bound on the quantum code distance of a quantum fractal liquid 
can be obtained by considering the following polynomial decompositions: 

7(a:)g(a;)^ = 5]a:7(xFCg)(7) (48) 
l{x)f{xY = Y,x^g{x)W<ifij) (49) 

ij 

where j{x) is an arbitrary polynomial and tensors C^j\'j) and clj\'y) have entries over Fp. Then one has 

d>mm{W{C{j)),W{D{j))} 

where W{C{j)) and W{D{j)) represent weights of tensors c'^j\j) and £)j-^'(7). 

While one needs to consider all 7(2;) to obtain a lower bound, tensors C^j\j) and dIj\j) with minimal weights 
often appear for 7(2:) = 1. One needs to consider the following decompositions: 

[-(x-2)]^ = ^:rM^(:r + 2)FC(f (50) 

ij 

[^(a; + 2)]^ = ^:rM-(x-2)Fi?(f. (51) 

ij 

By absorbing constant factors and noting symmetries between (a; — 2) and (a; + 2), it suffices to consider the following 
decomposition: 

{x~lY^^x\x + iyC^,f (52) 

ij 

^ (£) . 

where C^^j is a renormalized tensor over Fp. Noting that 

(X - lY =x'-l- .T^-l + ,C2 • _ . ^£-3 . . . 

{x + \f = x'^l-x'-^^,C2-x'-^^tC^-x'-^--- 

the following decomposition may be optimal 

(x ~ \Y = (x + 1)^ - 2 • fCi • x^-^ ~2-iC3- x^-^ ■■■ . (53) 

We performed numerical simulations for p < 7, and found that the above decomposition is indeed optimal. 
Assuming that the above claim holds for p > 7, one obtains 

- Z^i 2 ' 8 

where \n\ represents the smallest integer which is larger than (or equal to) n. One obtains a local quantum code with 

n-0(p3), d-0(/), fc(P)-0(p). 
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By simulating it with two-dimensional spins, one has 

n^p^logp, d^p^logp, k ^ plogp 

where the range of interaction terms is w logp which is logarithmically divergent. So, quantum fractal liquids may 
saturate the local code bound in three spatial dimension with 

fc~0(L), dr^O{L^), wr^O(\ogL). 
VII. CONCLUSIONS 

In this paper, we have presented general construction of classical and quantum fractal liquids and discussed their 
physical properties. We have found that correlation functions of classical fractal liquids have peculiar oscillatory 
behaviors which correspond to power-law decay with imaginary scaling dimensions. We have demonstrated that 
degenerate ground states of classical and quantum fractal liquids have discrete scale symmetries and correspond to 
limit cycles under real-space RG transformations. We have argued that quasi-particle excitations in classical and 
quantum fractal liquids violate charge conservation as they possess algebraic symmetries which are more general 
than gauge symmetries. We have derived a necessary and sufhcient condition for quantum fractal liquids to be free 
from string-like logical operators by using an analogy to cellular automata. We have demonstrated that topological 
entanglement entropies can detect quantum fractal order in principle. We hope that our construction and theory of 
fractal spin liquids will provide a stepping stone toward complete understanding and classification of quantum phases 
of matter. 

Quasi-particle excitations arising in quantum fractal liquids are all Abelian. Topological order arising in quantum 
fractal liquids is non-chiral. It is unclear to what extend our results may generalize to chiral topological phases and 
non- Abelian topological phases. Effective field theoretical descriptions of classical and quantum fractal liquids are 
currently not known. It may be interesting to analyze how classical and quantum fractal liquids behave under RG 
transformations in the language of matrix and tensor product state representations [55-n58j . Connection to limit 
cycles in non-relativistic field theory has not been established. Classification of quantum phases arising in translation 
symmetric stabilizer Hamiltonians is an important problem. In [43], we proved that topological phases arising in 
two-dimensional translation symmetric Hamiltonians are equivalent to copies of the Toric code. Yet, classification 
of three-dimensional topological phases remains as an open problem. An underlying difficulty in physically realizing 
the Sierpinski triangle model lies in the fact that the model has three-body interaction terms. Yet, one may simulate 
three-body interactions easily by using hopping particles as proposed in |59) . Note that fractal geometries appear in 
various biological systems, such as broccoli and sea shells in nature. So, we think that physical realization of fractal 
geometries is possible in principle, and is within the reach of current experimental techniques. 

Conventional classical spin liquids are systems that violate the third law of thermodynamics with extensive ther- 
modynamic entropy at zero temperature. Whether one can have topological spin liquids that violate the third law of 
thermodynamics is an interesting open problem. According to an analytical result obtained in quantum information 
community |52) . for commuting frustration- free quantum spin systems, one cannot have topological spin liquid that 
violate the third law of thermodynamics. Yet, the problem is open for frustrated systems. Indeed, frustration is the 
very origin that leads to all the interesting phenomena in both gapped and gapless spin liquids. 
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Appendix A: Tool-box for translation symmetric stabilizer Hamiltonian 



Translation symmetric stabilizer Hamiltonian is a particularly interesting class of exactly solvable quantum spin 
systems which may possess rich varieties of quantum phases. Systematic studies of translation symmetric stabilizer 
Hamiltonians were initiated in |34|, I43j where topological phases of two-dimensional stabilizer Hamiltonians are com- 
pletely classified and are shown to be equivalent to copies of the Toric code. It has been also pointed out that 
translation symmetric stabilizer Hamiltonians can be studied via commutative algebra [25) . 

In this appendix, we develop theoretical tools for analyzing physical properties of translation symmetric stabilizer 
Hamiltonians. For simplicity of discussion, we concentrate mainly on Hamiltonians with the following interaction 
terms: 




where a and /3 are arbitrary polynomials over Fp. We refer to the model as canonical model. 



Logical qubits: Consider a polynomial R over Fp which satisfies 

Ra = RP = 0. (A2) 
The number of logical qubits corresponds to the number of solutions to the above equation since R generates an 



identity operator: Z{Ra, R/3)'^ = I. We call such polynomials R the identity generators and call Eq. (A2) the 
self- consistency equation. The total number of logical qubits is 

fc = 2dimi? (A3) 

by counting contributions from Z-type and X-type stabilizer generators. Note there are p'' ground states in total. 
For instance, with a = 1 + and (3 — 1 + over F2, one needs to find solutions to i?(l + x^) = i?(l + x) = 0. 
Note 1 -I- x'^ = (1 -I- a;)(l -|- a; -I- x"^). For L = 2™, one finds that (1 -f a; + x"^)^ = 1, and {\ + x + x^) has an inverse 
under periodic boundary conditions. Then the self-consistency equation is further reduced to i?(l -|-a:) = 0. Note that 
(1 -I- x) is not reversible since (1 + x)^ = 0. The only possible non-trivial solution to the self-consistency equation is 

i? = (1 + x)^^^ = l + x + x'^ A 1- and one has dimi? = 1 and k = 2. 

The irreversibility of a and /? is important in finding identity generators R. Let us assume that exists. Then 
the self-consistency equation for R is given by Ra = which has i? = only as a solution. So, in order for A: ^ 0, a 
and P need to be irreversible. It turns out that irreversibility of a and /? is a necessary and sufficient condition for 
fc 7^ when L — p™. Recall that, for a polynomial / = c,iX% one has P = CiX^^ . So, for L — p™, one has 



One may notice the following: 



/ is irreversible <^ f{x = 1) = 0. (A5) 



Similarly, for higher-dimensional systems, one has 

/ is irreversible f{x = 1, y = 1, 2; = 1, • • • ) = 0. (A^6) 

Then, R = ^ ^ x'^y^ ■ ■ ■ is an identity generator since Ra — Rf3 — for irreversible a and /3. Therefore, one has 
the following criteria for logical qubits: 

A model has a logical qubit (fc 7^ 0). a and f3 are irreversible. 

Topological order: A necessary and sufficient condition for a canonical model to have topological order can 
be concisely described by algebraic relations between a and p. Without loss of generality, one may assume that 
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a and /3 have terms with positive finite powers oi x,y, z ■ ■ ■ only. (Otherwise, we shift the lattice positions). We 
also assume k ^ Q. For now, let us lift periodic boundary conditions and consider a and /3 on a semi-infinite lattice 
(a;, y, z, • • • = 0, • • • , go). Let us assume that a and /3 have a common factor 7: 

a = 7ao, /3 = 7^0 

on a semi-infinite lattice. Here, we consider the largest common factor 7 between a and /3. Note that ag, /3o and 7 
are finite polynomials with local supports since a and /3 are finite polynomials due to the locality of interaction terms. 

Now, let us return to a system under periodic boundary conditions {x^ = 1). Then, one has the following criteria 
for topological order (assuming fc 7^ 0): 

A canonical model has topological order. <^ The largest common factor 7 of a and /? is reversible. 

We emphasize that the reversibility of a common factor 7 is defined under periodic boundary conditions while a 
common factor 7 itself is derived on a semi-infinite lattice. 

It is easy to prove "=4>" of the above claim. When a and /3 have an irreversible common factor 7, one can show 
that the following operator is a logical operator: 

which leads to d ^ 0{1). To prove "<J=" of the claim, we assume that a canonical model does not have topological 
order and there exists a finite logical operator i^^'' = Z{a, b)^ with finite polynomials a and b with a(3 — ba = 0. 
Suppose that 7 is reversible. Then Z{aQ, /3o)^ is a stabilizer, and a/3o — bag = 0. As ao and /3q do not have common 
factor, one must have a — Sao and b — S^q for some finite polynomial 5. Then £^-^^ is a stabilizer, which leads to a 
contradiction. 



Local unitary transformation: Here, we discuss how a canonical model is transformed into the other via a 
sequence of two-qubit control gates. We then show that the Cubic code is unitarily equivalent to a certain second- 
order quantum fractal liquid. 

Two-qubit gates can be characterized by a two-qubit Pauli operator y = Vi (8" V2 . Consider an arbitrary two-qubit 
Pauli operator U — Ui (E) 1/2- A two-qubit gate generated by V transforms U as follows: 

Ui(g)U2 UiV^^®U2V^^ (A7) 

where 



U^V^ = {-If'ViUi, U2V2 - {-irV2U2 

For instance, with V — Xi (E) X2, one has the following transformations: 

Zi, Z2 1 ^ f Z1X2, X1Z2 
Xi, X2 J \ Xi, X2 

These two-qubit gates may be viewed as generalizations of control-Z operation. One may see that transformations in 



Eq. (A7) preserve commutation relations by direct calculations, and thus are indeed unitary transformations. 

Let us apply these two-qubit gates to a canonical model with Z{a,f3)'^ and X{P,a)'^. We think of applying a 
sequence of two-qubit gates, characterized by Xa (8 Zb, on neighboring sites in the x direction. Then, one has the 
following transformations: 



with a — >■ a and /? — )• /3 -I- xa. Note that these two-qubit gates can be applied simultaneously as they commute with 
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each other. By generalizing this transformation, the fohowing transformations are allowed: 

a-^a, P^fi + x'-y^z^a (A9) 

where are some finite integers. 

For the Cubic code (a = 1 + x + y + z, j3 = l + xy + yz-\-zx)^ we apply two-qubit gates (a, /3) — >■ (a, /3 + xq;), modular 
transformations (x, y, z) — ^ (x, yz^^ , z), shifting of lattice sites in z directions and two-qubit gates (a, /3) — > (a + /3, /3): 

1 + x + y + z \ / 1 + x + y + z \ f 1 + x + yz~^ + z 
1 + xy + yz + zx J \l + x + x'^+yzj \ 1 + x + x'^ + y 

y + {l + x)z + z^\ I' {1 + X + x2) + (1 + x)z + z2 
1 + X + x^ + y J ^ \ 1 + x + x^ + y 

This is a second-order quantum fractal liquid. 

Open boundary condition: One may also consider open boundary conditions by defining stabilizer generators 
nicely so that they commute with each other at boundaries. Without loss of generality, one can assume that a and /3 
have the following forms: 

a = 1 + cix + C2X^ + • • • J P = 1 + dix + d2X^ -I- • • • . 
For Z-type stabilizer generators for j — Q, - ■ ■ , L — 1, we impose the following constraints: 

a;-' = for j > L. 



In other words, Z-type stabilizers are truncated at the boundary on the right hand side (Fig. 11). X-type stabilizer 
generators are defined by taking the duals of Z-type stabilizer generators. They are truncated at the boundary on 
the left hand side. One may see that stabilizer generators commute with each other at the boundaries. A similar 
construction works for higher-dimensional systems too. We note that this construction of canonical models with 
open boundary conditions reproduces a model of the Toric code with boundaries, originally proposed by Bravyi and 
Kitaev where boundaries are classified into two types; rough and smooth boundaries [60^ . 

Higher-dimensional model: The /^-dimensional Toric code with m-dimensional and (Z? — m)-dimensional logical 
operators can be described by a model with oCm qubits per site with variables xi, • • • ,xo- We label sublattices by 
a set of m integers {Ai,A2, ■ ■ ■ ,Am) for Aj E Z/j. For Z-type stabilizer generators, we choose an arbitrary set of 
m+1 integers A* from Z^. If {Ai,A2, • • • , A^) U Aq = A* , we put 1 -I- xaq for a sublattice {Ai,A2, • • • , Am). As for 
X-type stabilizer generators, we consider a dual of {Ai, A2, • • • , Am) (a maximal set of integers which do not overlap 
with {Ai,A2, • • • , Am)) and choose an arbitrary set of D — to + 1 integers. 

For instance, with D = 3 and m = 1, three sublattices are labeled by (1),(2),(3). One may consider A* ~ 
(1, 2), (2, 3), (3, 1) and obtain Z-type stabilizers: 




As for X-type stabilizers, sublattices are labelled by duals (2, 3), (3, 1), (1, 2). Taking A* = (1, 2, 3), one has 




With D = 4 and TO = 2, sublattices are labelled by (1, 2), (1, 3), (1, 4), (2, 3), (2, 4), (3, 4), and for A* = 
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FIG. 11: A construction of stabilizer Hamiltonians with open boundary conditions. An example with a = 1 -\- x + and 
/? = 1 + a; is shown. For Z-type stabilizers, Pauli operators are truncated at the right edge, and for X-type stabilizers, Pauli 
operators are truncated at the left edge. 



(1,2,3), (1,2,4), (1,3,4), (2,3,4), Z-type stabilizers are 
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As for X-type stabilizers, sublattices are labelled by duals (3, 4), (2, 4), (2, 3), (1, 4), (1, 3), (1, 2), and one has 



/ \ 



l+Xi 



1+^2 
V 1 + ^3 / 
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l+Xl 
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V 1 + ^4 y 
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To construct fractal versions of these models in {D + 1) -dimensional systems, we add an extra variable xq, and replace 

Xj fj{xo)xj 

where fj{xo) is a polynomial with xq only. 

Lower bound on code distance: Below, we reduce a problem of finding lower bounds on code distances of 
quantum fractal liquids to a certain problem of decomposing polynomials f{x) and g{x). Consider a quantum fractal 
liquid for L — with reversible f{x) and g{x) with f(x = 1) = 1 and g{x = 1) = 1. Then, one has the following 
holds. 
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Consider all the tensors C,y satisfying the following equation: 

-/ix)g{xY = J2^'f(^y^l?^ fo^all L (AlO) 

ij 

and consider all the tensors D^^^ satisfying the following equation: 

-/{x)f{xY = ^x'g{xyD'^f, for all t (All) 



We denote minimal weights of tensors C^f and D-j^ as Wm^ and wjn'^ for 7(0;) ^ 0. We denote the minimum 

e dm is lower bounded by 

dm > Wm. (A12) 



of Wm^ and as Wm- Then the code distance dm is lower bounded by 



Below, we present the proof of the above claim. Consider a quantum fractal liquid with L x L x L sites with 
L = p™, and denote the ^-th layer of L x L x 1 sites as for ^ = 0, • • • , i — 1. We think of removing sites in a subset 
E = EqV} Ei\J ■ ■ ■ \J El-i with E being a set of site indices and Ei C Q^. We ask if £(^)-typc logical operators can 
be supported inside the remaining sites E or not. ^^^^-type logical operators can be represented as: 

^(-) (7) = ^ ( ^^J^^^ ) , 7(x) = CO + ci. + . . . + c,.,x^-^ 

where Cj = 0, • • - p — 1. Let us denote coefficients of x'y^ in 'y{x)f{x,y) as Sijo{'y)- 

-f{x)fix,y) = J2Sijoh)x'y^- 

ij 

Observe that a logical operator i^^^-y) can be supported inside Qo/{i,j, 0) if and only if a polynomial y{x)f{x, y) has 
a zero coefficient for x'^y^ . Then £(^^(7) is supported inside Qo/Eq if and only if S'ijo(7) = for all {i,j,0) € Eq. 
It is convenient to introduce polynomials Rijo which satisfy the following relations: 

R^^, = r(^^°) + rf^°)x + ■ ■ ■ + rE°^^-\ where SM = E r^^^c^. (A13) 

a=0 

For instance, one has i?ooo = 1, -Rioo = x and i?2oo = a;^ since 5ooo(7) = cq, <S'ioo(7) — ci and 5*200(7) = C2- With 
some speculation, one can prove that 

i?,,o = x'fixy (A14) 

where f{x) is the dual of f{x). Rijo may be viewed as a polynomial representation of coefficients S'ijo(7)- In order for 
£^^Hl) to be supported inside Qq/Eq, 7(0;) needs to satisfy a set of constraints that are characterized by polynomials 
inside S{Qo) = ( Rijo ■ («, j,0) e ^0 )• 

Next, let us find logical operators which have equivalent representations inside Qi/Ei. Note the following equiva- 
lence relation: 



^ ( 7(x)f(x, y))^^{ lix)9ix)fix, y)z)- 
So, £^^\j) can be supported inside Qi/Ei if and only if Siji{-)) = for all (i,j, 1) e Ei where Siji corresponds to 
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coefficients of x^y^z in z^{x)i{x^y)g{x). These constraints can be represented as 



L-l 



a=0 



In general, for polynomials 



one has 



where 



(A15) 



(A16) 



7(a:)f(x, ?/)g(a;, 2) = ^ S^-jtx'y z 

Consider the following spaces of constraints on ^{x) for Qi/Ef. 

SiQe)^{R^Ji■.{lJ,i)eEe). (A17) 
Then, all the logical operators can be supported inside E if 

s{Qo)ns{Qi)n---nS{QL-i)^<ii. 

In other words, if there is no common constraint for 7(0;) to satisfy among S{Qe), logical operators ^^'^■'(7) with 
arbitrary 7(0;) can be supported inside E. 

To analyze properties of E which satisfy Eq. ( A17), let us consider the case where S{Qe) have some common element 

(A18) 



for all t. Then, a tensor satisfying the following equation must exist 



-f{x) = ^x'f{xyg{xYcf^, foraU I 



Consider all the possible solutions c\f to Eq. (A18) and their weights (the number of non-zero entries), and denote 



the minimal weig ht as W!^n\l). Then£(^)(7) can be supported inside E if < Wrn (7)- Let us further define the 
minimum of Wm (7) for all 7(0;) 7^ as 

I^i^) =minW^if)(7). 

Then all the ^'^■^■'-type logical operators can be supported inside E if \E\ < w!^\ One can repeat a similar analysis 
for r'^^^-type logical operators, and ends up with the following equation 



(A19) 



By considering all the possible solutions to Eq. (A19) for all 7 7^ and denoting their minimal weight as wjn^\ one 
notices that all the r^'^^-type logical operators can be supported inside E ii E < Wr[^\ 

We denote the minimum of Wm ^ and as Wm- Then, all the logical operators can be supported inside E if 

l^-l < Wm- Then there is no logical operator whose weight is smaller than Wm [IS], and one has 



dm > Wm 



(A20) 
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where dm represents the code distance for L = p™. This leads to the proof of Eq. ( A12) 
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